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ABSTRACT 

We present an analysis of multi-epoch observations of the Type II-P supernova SN 2004et. 
New and archival optical spectra of SN 2004et are used to study the evolution of the Ha and 
[O I] 6300 A line profiles between days 259 and 646. Mid-infrared imaging with Michelle 
on Gemini North and with all three instruments of the Spitzer Space Telescope was carried 
out between 2004 to 2010, supplemented by archival Spitzer data. We include Spitzer 'warm' 
mission photometry at 3.6 and 4.5 fim obtained on days 1779, 1931 and 2151, along with 
ground-based and HST optical and near-infrared observations obtained between days 79 and 
1803. Multi- wavelength light curves are presented, as well as optical— infrared spectral en- 
ergy distributions (SEDs) for multiple epochs. Starting from about day 300, the optical light 
curves provide evidence for an increasing amount of circumstellar extinction attributable to 
newly formed dust, with the additional extinction reaching 0.8-1.5 magnitudes in the F-band 
by day 690. The overall SEDs were fitted with multiple blackbody components, in order to 
investigate the luminosity evolution of the supernova, and then with Monte Carlo radiative 
transfer models using smooth or clumpy dust distributions, in order to estimate how much 
new dust condensed in the ejecta. The luminosity evolution was consistent with the decay 
of 56 Co in the ejecta up until about day 6 90, after wh i ch an additional emission source is 
required, in agreement with the findings of Kotak et al. (2009). Clumped dust density distri- 



butions consisting of 20% amorphous carbons and 80% silicates by mass were able to match 
the observed optical and infrared SEDs, with dust masses that increased from 8 x 10 -5 M on 
day 300 to 1.5 x 10 -3 M on day 690, still significantly lower than the values needed for core 
collapse supernovae to make a significant contribution to the dust enrichment of galaxies. 
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1 INTRODUCTION 



Theoretical studies b y |Kozasa et al.| (|1989|>; |Todini & Fe rrara 
( |200!) ; |Nozawa et al.|p003| >; |Bianchi & Schneider| ( |2007| ) of dust 
formation in the ejecta of core-collapse supernovae (CCSNe) pro- 
duced by massive stars suggested that they could each produce up 
to 0.1-1.0 solar masses of dust and thereby act as major sources 
of dust in galaxies. The discovery at submillimetre wavelengths 
of large quantities of dust in very young high redshift galaxies, 
(e.g. |Bertoldi et al.|2003| |Maiolino et al.|2004| |Dwek et al.|2007) 
seemed to reinforce the case for massive star CCSNe as major dust 
contributors. However studies at infrared wavelengths of Type II 



CCSNe in the nearby Universe have to date not confirmed these 
dust production predictions. SN 1987 A in the LMC condensed 
of dust in its ejecta |Wooden et al.|l993* 



less than 10 - " 3 M 







Er- 



|colano et al.|2007) , while SN 2003gd in M 74 produced less than 
a few x 10 M of new dust by day 499 ( Sugerman et al.| 2006 



Mei kle et al.|2007| >. For a number of other recent Type II CCSNe 
estimates for the quantiti es of new dust formed in their ejecta have 
not exceeded 10" 3 M ( |Andrews et al.|2010l|201l[|Szalai et al. 



|20lT|pelkreetal.|2011| ). 

Kotak et al. (2009) studied the multi-wavelength evolution of 
the Type II-P SN 2004et, the eighth supernova to be discovered in 
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the past century in the nearby spiral galaxy NGC 6946. We report 
here our own multi-wavelength observations and analysis of the 
dust production by this supernova. Following |Li et al.| ( |2005| ), we 
adopt an explosion date of 22.0 September 2004 (JD 2,453,270.5) 
as the day epoch to which our observation dates are referenced. 
Section[2]provides a brief overview of previously published optical 
and infrared (IR) observations of SN 2004et. Section [3] describes 
the Ha and [O I] 6300 A line profiles obtained from a sequence of 
optical spectra of SN 2004et, while Section [4] describes our mid- 
IR observations. The resulting mid-IR fluxes are presented in Sec- 
tion [42] together with a discussion of the mid-IR light curve evo- 
lution. Complementary photometric data, obtained by us at optical 
and near-infrared (NIR) wavelengths using the HST and Gemini 
North telescopes, are presented in Section [5] The spectral energy 
distributions (SEDs) of SN 2004et at different epochs are discussed 
in Section [6] This section includes blackbody fitting to the SEDs, 
allowing limits to be placed on the physical properties of the SN 
at the various epochs. Using Monte Carlo radiative transfer mod- 
els, we estimate ejecta dust masses in Section [7] Finally, the re- 
sults from our multi-epoch analysis of SN 2004et are summarised 
in Section [8] along with the implications for dust production by 
Type II SNe. 



2 PREVIOUS STUDIES OF SN 2004ET 



2.1 Optical observations 



|Z witter et al.| |2004 ) reported that an echelle spectrum taken on 
28 September 2004 showed a relatively featureless spectrum with 
very broad, low-contrast Ha emission, suggesting SN 2004et to be 
a Type II event. This was confirmed by a low-resolution spectrum 
taken on 1 October 2004 ( Filippenko et al. 2004 ), which showed the 
P-Cygni profile of Ha to be dominated by the emission component, 
but that the other hydrogen B aimer lines had more typical P-Cygni 
profiles. The continuum was reported as quite blue but dropping off 
steeply shortwards of 4000 A. 

The V, R and /-band light curves showed SN 2004et to be a 
Type II-P (plateau) supernova ( |Li et al.|20 05 ) but it exhibited some 
differences when compared with the typical Type II-P SN 1999em. 
Evolution was slower in the U and B bands, consistent with the (U— 
B) and (B - V) colour evolution, leading Li et al. ( 2005) to conclude 
that SN 2004et seemed to evolve more slowly than SN 1999em, 
especially in the violet part of the spectrum. 

Extensive photometric and spectroscopic monitoring in the 
optical was carried out by Sa hu et al.| ([2006) from days 8 to 541 
after explosion. They confirmed the supernova to be of Type II- 
P based on the plateau observed in the VRI bands which lasted 
for ~ 110 days after the explosion. From their light curve analy- 
sis, they determined that the SN was caught at a very early stage 
soon after the shock breakout, reaching a maximum B-band mag- 
nitude ~ 10 days after explosion. They also found that the decline 
rate of the light curve during the early nebular phase (~ 180-310 
days) was similar to the radioactive decay rate of 56 Co, indicating 
that 7-ray trapping was efficient during this time, and estimated 
that 0.06 ± 0.02 M of 56 Ni was synthesised during the explo- 
sion of SN 2004et, in agreement with the subsequent estimates of 
|Misra et al!| ( |2007| ) and jMaguire et al.| ( |20T0] >. Based on the plateau 
luminosity and duration and the mid-plateau expansion velocity, 
|Sahu et al.| ( [2006%|Misra et al.|{2007) and |Maguire eTaL] pOTO} 
have estimated explosion energies for SN 2004et in the range (0.88- 
1.20) xlO 51 erg. 



2.2 Distance and reddening 

There are many estimates in the literature for the distance to the 
host galaxy NGC 6946. Sahu et al. (2006 ) summarised a few of 
these and included the result from their own analysis of SN 2004et 
using the standard candle method of Nuge nt et al.| (2006 ) for 
Type II-P supernovae, deriving an average distance of 5.6 Mpc. 
Herrmann et al.| ( |200 8) estimated a distance of 6.1 ±0.6 Mpc to 
NGC 6946 using the planetary nebula luminosity function. For con- 
sistency with our previous studies of SN 2002h h ([Barlow et al.| 
|2005l |Welch et al.||2007] ) and SN 2008S ( [Wesson et al.||2010) , 
we adopt the distance of 5.9±0.4 Mpc to NGC 6946 estimated by 
|Karachentsev et al.|(|2000 1 



From |Schlegel et ah 



from the brightest stars method. 



(1998 ), the foreground Galactic redden- 
ing towards SN 2004et is estimated to be E(B-V) = 0.34 mag. 
Z witt er et al.| |2004 ) used the equivalent width of Na I D 2 lines 
from their high-resolution echelle spectra to estimate a total red- 
dening (Galactic + host) towards the SN of E(B-V) = 0.41 mag. 
|Sahu et ah] ( |2006| > used similar analysis with their lower-resolution 
spectra to obtain a comparable value of E(B - V) = 0.43 mag. We 
adopt a total reddening to SN 2004et of E(B- V) = 0.41 mag. 



2.3 Evidence for dust production by SN 2004et 

From their spectroscopic observations of the temporal evolution of 
the Ha and [O I] 6300, 6363 A line profiles from days 277 to 465, 
Sahu et al. ( 2006 ) noted a blueshifting after day 300 of the emission 
peak of both features and a flattening of the Ha emission peak. 
This, together with a steepening of the light curve after day ~ 320, 
they interpreted as indications of early dust formation in the ejecta 
of SN 2004et. [Misra et al.| ( [2007] > also found that the rate of decline 
in the optical light curves accelerated between ~ 320-386 days. 

Fabbri, Sugerman & Barlow (2005) reported the day 64 detec- 
tion of SN 2004et in all four Spitzer Infrared Array Camera (IRAC) 
bands, from 3.6 to 8.0 pm, in Spitzer Infrared Nearby Galaxy Sur- 
vey (SINGS) Legacy program archival images of NGC 6946 (see 
Section |47TJ. Kotak et al. (2009 ) presented archival and their own 
Spitzer mid-IR observations of SN 2004et obtained up to day 1406, 
together with late-time optical spectra. They concluded that be- 
tween days 300-795, the spectral energy distribution was com- 
prised of three components - hot, warm and cold - each respec- 
tively due to emission from: optically thick gas; newly-formed, 
ejecta-condensed dust; and an IR echo from the interstellar medium 
of the host galaxy. They estimated that the mass of dust formed 
in the ejecta grew to a few times 10 -4 M©, located in co-moving 
clumps of fixed size. From their Spitzer Infrared Spectrometer 
(IRS) spectra, they reported the first spectroscopic evidence of sil- 
icate dust formed in the ejecta of a supernova, supported by the 
detection of strong but declining molecular SiO emission in the 8- 
//m region. They interpreted the appearance of broad, box- shaped 
optical emission line profiles about 2 years post-explosion as due to 
the impact of the ejecta on the circumstellar medium of the progen- 
itor star, resulting in the formation of a cool, dense shell to which 
they attributed responsibility for a later rise in the mid-IR emission 
fro m SN 2004et. 

Mag uire et al.| < |2010) > reported optical and NIR photomet- 
ric and spectroscopic observations of SN 2004et carried out from 
just after explosion to +500 days. Their NIR spectrum at day 
306 showed a clear detection of the first overtone band of CO at 
^2.3 /xm, which they interpreted as being a signature of dust for- 
mation. By analysing the optical light curves in the early nebular 
phase, they found that the BVR decline rates between days ~ 136 
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and 296 were consistent with those expected from light curves pow- 
ered by the radioactive decay of 56 Co (assuming complete 7-ray 
trapping). However, between days ~ 296 and 414 they found that 
these decline rates had steepened. They also noted a significant 
blueshift in the peak of the Ha emission line from days ~ 300 to 
464, in agreement with the results of Sahu et al. (2006 ). These re- 
sults were interpreted as signatures of dust formation occurring post 
300 days. Their very late time (> 1000 days) photometry showed 
a flattening of the optical and NIR light curves, which they mainly 
attributed to the interaction of the SN ejecta with the circumstellar 
medium (CSM), following the work of |Kotak et aE] ( [2009l >. 



3 OPTICAL SPECTROSCOPY OF SN 2004ET 

As part of the SEEDS program, Gemini North spectroscopic obser- 
vations of SN 2004et were obtained on 5 August and 3 1 October 
2005, corresponding to 317 and 404 days after explosion. A log of 
these and the other Ha-region observations that we analyse here 
can by found in Table [T] We utilised GMOS-N equipped with a 
0"75 slit width and the B600-G5303 grating in long-slit mode with 
a position angle of 296°. Three spectra were obtained during each 
epoch with identical exposure times of 900 s. The central wave- 
lengths of the images were 5950, 5970, and 5990 A, respectively, 
to allow for gap removal of the combined spectra. A 2 x 2 binning 
of the CCD pixels in the low gain setting was employed. All spectra 
were taken with adjacent GMOS baseline calibration flat exposures 
to correct for sensitivity gradients across the CCD, and CuAr spec- 
tra were utilised for the initial calibration of the dispersion solution. 

Our GMOS-N spectra were reduced using the Gemini IRAF 
package. Pipeline processed calibration images were obtained from 
the Gemini Science Archive. The spectra were trimmed and then 
overscan, bias and flat-field corrected using the task gs reduce. 
Wavelength calibration solutions were determined from the CuAr 
lamp spectra using gswavelength, and the solution was applied 
to the SN 2004et spectra via gstransf orm. Object spectra were 
extracted using gs extract. The observations were not flux cal- 
ibrated since our primary goal was to monitor the evolution of the 
line profiles. 

|Sahu et aL] < |2006] > presented photometric and spectroscopic 
data for SN 2004et from approximately 8 to 541 days after the 
explosion. Their results showed a shift to the blue of the central 
peak of both Ha and [OI] 6300,6363 A at late times. They con- 
cluded that this was indicative of new dust formed in the ejecta 
of SN 2004et. However, comparisons between the Ha profiles in 
our GMOS-N spectra and those at similar epochs in the archived 
calibrated spectra of Sahu et al. (2006), obtained from the Online 
Supernova Spectrum Archive (SUSPECT^ showed some incon- 
sistencies, with their spectra showing clear blueshifting of the Ha 
emission line profile between days 314 and 391, while our own 
spectra showed little change between days 317 and 404. In order 
to help resolve this discrepancy, the Sahu et al. ( 2006) authors gen- 
erously provided their raw data for SN 2004et for several of the 
epochs presented in their paper. This allowed us to reduce the re- 
spective sets of data in the same way. A description of the observa- 
tional setup used by Sahu et al. ( 2006) can be found in their paper. 

We also obtained two archival spectroscopic observations of 
SN 2004et, taken on days 336 and 646. The earlier epoch cor- 
responds to an observation on 29 August 2005 taken with the 



http://bruford.nhn.ou.edu/~suspect 



Low Resolution Spectrograph on the 3.58 -m Telecopio Nazionale 
Galileo (TNG) on La Palma. The latter epoch corresponded to an 
observation on 30 June 2006 with the Faint Object Camera and 
Spectrograph (FOCAS) on the 8.2-m Subaru telescope on Mauna 
Kea in Hawaii. The reduction procedure matched that carried out 
by us for the GMOS-N and |Sahu et allpOOoT ) spectra. 

Since our goal was to self-consistently align a number of spec- 
tra taken with four different instrument/telescope setups, we did not 
rely solely on the initial wavelength calibrations. During the extrac- 
tion of each SN spectrum, a sky spectrum was also extracted that 
was ultimately subtracted from the SN spectrum. Since the sky and 
SN spectra possessed identical wavelength calibrations, we derived 
corrections to the initial wavelength calibrations of each of the SN 
spectra using the strong and narrow [O I] 5577 A and 6300 A sky 
emission lines. 

The temporal evolution of the Ha profile between days 259 
and 646, following our re-analysis of the combined set of spectra, 
is depicted in the left-hand panel of Figure [T] The spectra are dis- 
played in velocity space, with the continua subtracted and the peaks 
of the emission line normalised to approximately unity. The inset 
shows a close-up view of the earliest and latest Ha profiles, ob- 
tained just over a year apart at days 259 and 646, to highlight the 
overall blueshifting of the profile during this time. The dashed lines 
in the inset indicate the line peak centres measured on days 259 and 
646, with the blueshifting of the latter profile providing strong evi- 
dence for the formation of dust in the ejecta during the intervening 
period, as described by |Lucy et ah] ( |1989| > for SN 1987A. 

The Ha line profiles presented by | Sahu et al.| |2006) showed a 
significant blueshifting with time of the emission peak, amounting 
to a few hundred kms -1 , with the largest shift appearing to oc- 
cur between days 277 and 314. |Sahu et al.|p0 06) did not quantify 
the shifts in their Ha profiles, beyond stating that a blueshift in the 
emission peak was clearly seen after day 300 in their da y 277-465 
Ha (and [O I] 6300,6363 A) profi les. |Kotak et al.|(|2009| > confirmed 
this from their own analysis of the Sahu et al. (2006) spectra, report- 
ing a shift of —400 kms -1 in the whole Ha profile between days 
301 and 314, but little sign of a progressive blueshifting during 

{2006} 



the subsequent day 314-465 period. Using the Sahu et al. 



spectra, along with additional spectra, Maguire et al. (2010) found 
that the peak of the Ha emission line was at +280 ± 50 kms -1 
between days 163 and 300, but from days 314 to 464 showed a 
constant blueshift, to — 1 37 km s - 1 . 

We have used our recalibrated spectral dataset to measure 
the wavelength of peak Ha emission in each of the profiles plot- 
ted in Figure [T] The evolution of the peak wavelength is shown 
in Figure [2] The line peaks show an overall blueshifting between 
days 259 and 646 of ~ 13.2 A, corresponding to a velocity shift of 
~ —600 kms -1 . Between days 259 and 314 we measure an initial 
blueshift of ~ 3 A (« 140 kms -1 ), not as large as the 400 kms -1 
shift between days 301 and 314 measured by |Kotak et al.| j2009) 
from the |Sahu et aLlp006] > spectra in the SUSPECT archive. While 
we find little change between days 314/317 and day 39 1 in the mea- 
sured emission line peaks in the recalibrated spectral dataset (Fig- 
ure [2]), between days 404 and 428 we measure a blueshifting of 
the emission peak by ^4 A (185 km s -1 ). followed by a further 
blueshifting of ~ 4.5 A (205 km s -1 ) between the spectra obtained 
at days 428 and 646. We therefore find that the majority of the 
blueshifting occurred after day 391. 

Whilst there is no strong evidence for a developing line asym- 
metry in the earlier epoch Ha profiles (Figure [TJ, a diminution of 
the red wing can be discerned by days 428 and 646. By day 646 
a significant blueshifting (~ 4.5 A) of the whole profile since day 
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Table 1. Log of optical spectroscopic observations of SN 2004et from the SEEDS collaboration (GMOS-N, PI: Clayton), from 
the authors Sahu et al. ( 2006 1 and from online archives (TNG and Subaru). 



Date 


Age 


Telescope/ 


Wavelength 


Exp. Program ID 


Principal 




[days] 


instrument 


range [A] 


time 


Investigator 



2005-06-07 
2005-08-01 
2005-08-05 
2005-08-29 
2005-10-17 
2005-10-31 

2005- 11-23 

2006- 06-30 



259 HCT HFOSC a 

314 HCT HFOSC a 

317 Gemini GMOS-N 

336 TNG LRS b 

391 HCT HFOSC a 

404 Gemini GMOS-N 

428 HCT HFOSC a 

646 Subaru FOCUS c 



3500-7000; 5200-9200 1 x 900 s 

3500-7000; 5200-9200 1 x 900 s 

3500-10000 3 x 900 s 

3890-8000 lx 1800 s 

3500-7000; 5200-9200 1 x 900 s 

3500-10000 3 x 900 s 

3500-7000; 5200-9200 1 x 900 s 

4670-8970 2 x 900 s 



GN-2005B-Q-54 
AOT12 CAT-G109 

GN-2005B-Q-54 

S06A-152 



Sahu et al. 


2006 




2006 


G. Clayton 
E. de la Rosa 


Sahu et al. 


2006 


G. Clayton 


Sahu et al. 


2006 



a Raw data, including calibration frames and flux standards, provided by the authors i Sahu et al. 2006 . 
b Raw and calibration data were downloaded from the online TNG archive at http://ia2.oats.inaf.it/. 

c Raw and calibration data downloaded from the online Subaru Mitaka Okayama Kiso Archive (SMOKA) at http://smoka.nao.ac.jp/index.jsp. 
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Figure 1. The Ho; and [O I] profile evolution. The left hand panel shows the Hen spectra from days 259 to 646, with the normalised continua subtracted 
to a zero-level and the peaks of the emission line normalised to unity. The dashed vertical line in the main plot corresponds to the rest wavelength of Ha 
(NGC 6946 has a radial velocity of +40 km s _1 ). The inset shows a comparison of the earliest and latest Hen spectra, obtained just over a year apart, to 
highlight the blueshifting (~ 600 km s _1 ) during this interval. The dashed lines in the inset plot indicate the measured line centres. The right hand panel shows 
the evolution of the [O I] 6300 A profile from days 259 to 428. The dashed line in the main plot indicates the rest velocity of the 6300 A l ine. The emission 
feature in the red wing corresponds to the [O I] 6363 A line. The spectra from days 259, 314, 391 and 428 are from |Sahu et al.||2006) (raw data kindly 
provided by the authors and calibrated by us). The spectra from days 317 and 404 are SEEDS GMOS-N observations, and those from days 336 and 646 are 
from archival TNG-LRS and Subaru-FOCAS spectra respectively. 



428 is evident, with multiple peaks evident at the centre of the pro- 
file, although there are earlier inflections in the blue wing of the Ha 
profiles. Two broad features either side of the main emissio n line 
(at ~ +7000 and -6000 krns" 1 ) have appeared by day 646. |Kotak 
|et al.| ( |2009| > presented three later optical spectra of SN 2004et, ob- 
tained with the Keck telescopes on days 823, 933 and 1 146, which 
show similarities to, and a development of, the Subaru day 646 
Ha profile. They described the Ha profile from the late-time Keck 
spectra as having a steep- sided, box-like component, with a half 
width at zero intensity (HWZI) of 8500 km s _1 and noted that the 
characteristic ejecta profile of Ha seen in the |Sahu et al.| <j2006 ) 
spectra may still be present at days 823 and 933, with a HWZI of 
~ 2000 km s _1 , but that its presence at day 1 146 was less certain. 

The 6300 and 6363 A lines of [O I] can be seen in emission 
on the blue side of the Ha profiles plotted in Figure [T] In the right- 
hand panel of Figure [I] we have plotted the [O I] 6300- A velocity 



profiles for days 259 through to 428 (the day 646 profile had too 
low a signal to noise). The inset compares the day 259 and day 428 
[O I] profiles, showing a clear diminution of the red wing of the day 
428 profile relative to day 259. 



4 MID-INFRARED OBSERVATIONS 

4.1 Gemini-Michelle and Spitzer photometry 

NGC 6946, the host galaxy of SN 2004et, was observed with 
Spitzer by the SINGS Legacy program ( Ken nicutt et al.|2003] > be- 
tween June and November 2004, such that the region of the SN was 
serendipitously imaged pre- and post-explosion. Figure [3] shows 
the position of the SN in relation to its host galaxy NGC 6946 in 
the pre- and post-explosion SINGS IRAC images, together with a 
closer view of the SN region in the different IRAC wavebands. 
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Figure 3. The position of SN 2004et in SINGS Legacy IRAC images of NGC 6946. Panels (a) and (b) show the whole galaxy at 3.6-/im pre- (day -104) and 
post-explosion (day 64) respectively. The square is centred on the SN coordinates, with the SN clearly evident at day 64. Panels (c), (d) and (e) zoom in on the 
square region of the SN field at 4.5, 5.8 and 8.0 /im respectively, at day 64. 
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Figure 2. The evolution with time of the measured wavelength of the Hen 
emission peak of SN 2004et. 



SN 2004et was clearly detected in the SINGS IRAC image at 
day 64 and was subsequently monitored with IRAC, MIPS and the 
IRS Peak-Up Imaging (PUI) module via our Spitzer GO programs 
during Cycles 2, 3, 4, 6 and 7 (Cycle 6 was the beginning of the 
Spitzer 'warm' mission, in which the observatory operates using 



only the 3.6 and 4.5 /im IRAC channels). Archival Spitzer data for 
the SN from Cycles 2 and 3 were also downloaded to provide a 
more complete time- sample of the SN's mid-IR evolution. In sum- 
mary, post-explosion mid-IR observations of SN 2004et were taken 
between November 2004 and August 2010, corresponding to an age 
range of 64 to 2151 days. With the Spitzer Space Telescope, there 
are 12 epochs of observations in each of the four IRAC wavebands, 
with a further three epochs of IRAC observations during the 'warm' 
mission in just the 3.6- and 4.5-/xm channels. There are 8 epochs 
of IRS PUI observations at 16 /mi, and 9 epochs of MIPS 24-/xm 
observations. Of the 33 individual Spitzer observations listed in Ta- 
ble [2] the first 3 were obtained by the SINGS Legacy program, 5 
were obtained by programs led by Pis Meikle and Kotak and 25 
were obtained by our SEEDS program. 

In addition to the Spitzer mid-IR observations, broad-band N' 
photometry at 11.2 fim was obtained with Michelle on Gemini- 
North during 2005, 2006, 2007 and 2008, consisting of six obser- 
vations at the 4 epochs. A time-ordered list of all the Spitzer and 
Michelle mid-IR observations of SN 2004et is provided in Table [2] 
The listed exposure time is the total time spent on- source. The fi- 
nal column provides a key to the list of observing programs from 
which these data were obtained; the program number and principal 
investigator are detailed in a footnote to the table. Data retrieved 
from the archive are marked as such. Appendix | A 1 1 describes how 
the mid-IR data were processed. 
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Table 2. Summary of the mid-infrared imaging observations of SN 2004et with the Spitzer Space Telescope and Gemini North. 
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2007-08.-02 


1044 


IRS-PUI 


16.0 




1.0 x 1.2 


1.2 


283.1 


L91 


2007-08-12 


1054 


IRAC 


3.6/4.5/5.8/c 


;.o 


5.2x5.2 


0.75 


321.6 


[9] 


2007-08-27 


1069 


MIPS 


24.0 




5.4 x 5.4 


0.75 


420.0 


[9] 


2007-12-09 


1173 


IRS-PUI 


16.0 




1.0 x 1.2 


1.2 


283.1 


[9] 


2007-12-27 


1191 


IRAC 


3.6/4.5/5.8/c 


;.o 


5.2x5.2 


0.75 


321.6 


[9] 


2008-01-07 


1202 


MIPS 


24.0 




5.4 x 5.4 


0.75 


420.0 


[9] 


2008-01-17 


1212 


IRS-PUI 


16.0 




1.0 x 1.2 


1.2 


1258.4 


[10] 


2008-06-21 


1368 


Michelle 


11.2 (NO 




0.5 x 0.4 


0.1 


1928.6 


[11] 


2008-07-09 


1386 


Michelle 


11.2 (NO 




0.5 x 0.4 


0.1 


2257.9 


[12] 


2008-07-18 


1395 


IRAC 


3.6/4.5/5.8/c 


;.o 


5.2x5.2 


0.75 


321.6 


[9] 


2008-07-29 


1406 


MIPS 


24.0 




5.4 x 5.4 


0.75 


420.0 


[9] 


2009-08-06 


1779 


IRAC 


3.6/4.5 




5.2x5.2 


0.75 


1161.6 


[13] 


2010-01-05 


1931 


IRAC 


3.6/4.5 




5.2x5.2 


0.75 


1161.6 


[13] 


2010-08-13 


2151 


IRAC 


3.6/4.5 




5.2x5.2 


0.75 


1161.6 


[14] 



II] 



Spitzer Cycle 1 SINGS Legacy program 00159, PI: Kennicutt. 



1 Archival data, Spitzer Cycle 2 GO program 20256, PI: Merkle 



^ This paper, Spitzer Cycle 2 GO program 20320, PI: Sugerman. 
[4] This paper, Gemini semester 05A program GN-2005A-Q-20, PI: Barlow. 
^ This paper, Gemini semester 06A program GN-2006A-Q-1, PI: Barlow. 
Archival data, Spitzer Cycle 3 GO program 30292, PI: Merkle. 



16] 



1 This paper, Spitzer Cycle 3 GO program 30494, PI: Sugerman. 

' This paper, Gemini semester 07A program GN-2007A-Q-5, PI: Barlow. 

' This paper, Spitzer Cycle 4 GO program 40010, PI: Meixner. 

)] Archival data, Spitzer Cycle 4 GO program 40619, PI: Kotak. 

L] This paper, Gemini semester 07B program GN-2007B-Q-4, PI: Barlow. 

' ] This paper, Gemini semester 08B program GN-2008B-Q-44, PI: Barlow. 

^ This paper, Spitzer Cycle 6 GO program 60071, PI: Andrews. Cycle 6 took place during the Spitzer 'warm' mission following com- 
pletion of the cryogenic mission. Only IRAC 3.6 and 4.5 /j,m channels were available, with expected sensitivity unchanged from perfor- 
mance in the cryogenic mission. 

^ This paper, Spitzer Cycle 7 GO program 70008, PI: Andrews, Spitzer 'warm' mission. 



4.2 Evolution of the mid-IR emission 



Table [3] lists the complete set of Spitzer and Gemini mid- 
infrared flux densities and associated uncertainties/upper-limits of 
SN 2004et from days 64 to 2151 as determined from the PSF-fitting 
techniques described by |Fabbri| ( [2011| ). Spitzer IRAC, MIPS and 
IRS-PUI data have had pre-explosion flux levels subtracted. Mea- 
sured (IRAC and MIPS) and estimated (PUI) pre-explosion flux 



densities at the position of the SN are summarised in the last row 
of the table. 

The IRAC, PUI 16-/xm and MIPS 24-/xm light curves (with 
pre-explosion levels subtracted) are shown in Figure [4] For clarity, 
the different light curves have been arbitrarily shifted by the fac- 
tors shown. Upper limits to the flux densities at 3.6 /am and 5.8 /am 
are indicated by the downward pointing arrows. Each waveband 
demonstrates the decline in brightness from the earliest epochs to 
around day 800 when the SN has faded or is fading to its faintest 
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Table 3. Mid-infrared photometry of SN 2004et. The shaded rows indicate those fluxes across different wavebands considered as the same epoch. 
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UT date 



Age 
[days] 



Flux density 1 / upper limits 2 [/dy] 



IRAC 3 



Michelle N' IRS-PUI 3 MIPS 3 



2004-06-10 


-104 


83.8 ± 10.2 


45.4 ± 8.3 


178 ± 25 


412 ±51 








2004-07-09 


-75 














376 ± 25 


2004-11-25 


64 


17490 ± 532 


13038 ± 745 


10046 ± 179 


6020 ± 91 








2005-07-13 


294 












930 ± 33 




2005-07-19 


300 


726 ± 61 


3151 ±97 


1291 ±218 


2162 ± 154 








2005-07-30 


311 










1700 ± 200 






2005-08-03 


315 














832 ± 75 


2005-09-17 


360 


430 ±71 


1728 ± 88 


935 ± 285 


1731 ± 156 








2005-09-24 


367 














735 ± 74 


2005-11-02 


406 


315 ±21 


1045 ± 21 


707 ± 79 


1500 ± 77 








2005-12-22 


456 












890 ± 32 




2005-12-30 


464 


174 ± 67 


656 ± 77 


606 ± 293 


952 ± 162 








2006-01-11 


476 














686 ± 89 


2006-05-12/14 597/599 










<650 






2006-08-04 


681 












670 ± 34 




2006-08-13 


690 


<49 


45.9 ± 15.3 


115 ±38 


342 ± 114 








2006-09-01 


709 














663 ± 77 


2006-09-10 


718 












562 ± 30 




2006-12-29 


828 


<45 


29.4 ± 9.8 


< 87 


145 ± 48 








2007-01-21 


851 














513 ±54 


2007-01-27 


857 












446 ±31 




2007-06-26 


1007 












324 ± 29 




2007-07-03 


1015 


<23 


14.1 ±4.7 


< 83 


113 ±37 








2007-07-09 


1020 










<400 






2007-07-10 


1021 














644 ± 48 


2007-08-02 


1044 












366 ± 35 




2007-08-12 


1054 


< 34 


15.1 ±5.0 


<78 


115 ±38 








2007-08-27 


1069 














610 ±49 


2007-12-09 


1173 












1023 ± 33 




2007-12-27 


1191 


52.5 ± 17.5 


189 ± 10 


345 ± 75 


458 ± 73 








2008-01-07 


1202 














1276 ± 42 


2008-01-17 


1212 












1055 ±31 




2008-06-21 


1368 










1036 ±212 






2008-07-09 


1386 










1016 ± 224 






2008-07-18 


1395 


76.3 ± 22.6 


258 ± 11 


465 ± 49 


578 ± 82 








2008-07-29 


1406 














1563 ± 54 


2009-08-06 4 


1779 


17.3 ±7.6 


113 ±8 












2010-01-05 4 


1931 


14.4 ± 6.9 


90.3 ± 16.8 












2010-08-13 4 


2151 


7.2 ± 3.0 


76.7 ± 5.3 












Pre-explosion 




83.8 ± 10.2 


45.4 ± 8.3 


178 ± 25 


412 ±51 




221 ± 22 


376 ± 25 



All post-explosion flux densities were measured with PSF-fitted photometry (using IRAF daophot). Pre-explosion IRAC and MIPS fluxes were measured in an 
aperture of radius 5" with sky annuli at inner radius l'.' 5 and outer radius 10 ;/ respectively, using a 2-cr clipped-mean sky algorithm (using IRAF phot). 

2 Flux upper limits for the non-detections in the Gemini-Michelle data are 3-cr values based on the standard deviation of the background in the region of the SN position 
scaled to a diffraction-limited size aperture. 

3 Spitzer IRAC, MIPS and IRS-PUI data have had pre-explosion flux levels subtracted: IRAC and MIPS by use of difference imaging techniques to subtract pre- 
explosion SINGS images, and IRS-PUI by estimating the pre-explosion level from blackbody fits to the SED at day 1015 (see text). Measured (IRAC and MIPS) and 
estimated (IRS-PUI) pre-explosion fluxes at the position of the SN are summarised in the last row of the table. 

4 Observations at days 1779, 1931 and 2151 were obtained during the Spitzer 'warm' mission, where only IRAC 3.6 and 4.5 /xm channels were available. 



levels. At 3.6 /xm and 5.8 /mi, the upper limits measured from the 
difference images indicate that the SN faded to below background 
levels for ~ 200 days. The distinctive rise in brightness after this 
time (> 1000 days) is evident in all wavebands. The latest Spitzer 
data at days 1779, 1931 and 2151 were obtained during the post- 
cryogenic phase of the mission, where only the shortest wavelength 
IRAC channels at 3.6 and 4.5 /xm were available. They show that 
sometime between days 1395 and 1779, the mid-IR brightness of 
the SN at 3.6 and 4.5 /im began to decline again, continuing with a 
slower decline to day 2151. The 3.6 /xm flux at day 2151 is about 



9 % higher than the pre-explosion level, whilst the 4.5-/mi flux at 
the latest epoch is a factor of 2.7 brighter than the pre-explosion 
level. Whilst the 3.6 /mi flux densities at days 1779 to 2151 are 
lower than the upper limits between days 690 and 1015, the on- 
source integration times for the later observations were between 
factors of 3 and 80 longer than those for the earlier observations. 
The corresponding increase in signal-to-noise for the IRAC 'warm' 
images, together with difference imaging techniques, allowed the 
SN to be reliably detected at deeper levels than for previous obser- 
vations with shorter exposure times when the SN was also faint. 
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Figure 4. Spitzer mid-infrared light curves of SN 2004et: 3.6, 4.5, 5.8 and 8.0 pm (IRAC), 16 pm (IRS Peak-Up Imaging) and 24 pm (MIPS). The IRAC and 
MIPS flux densities are from PSF-fitted photometry carried out on difference images which use the pre-explosion image (day -104 for IRAC and day -75 for 
MIPS) as the reference image. The 16 pm flux densities are from PSF-fitted photometry of the original images from which an estimated pre-explosion flux 
was subtracted (see text for details). For non-detections, upper limits to the flux densities are indicated by the downward-pointing arrows. For clarity the light 
curves have been shifted vertically by the factors indicated. 



A selection of pre- and post-explosion IRAC images at 
4.5 pm, 5.8 pm and 8.0 pm are shown in Figure [5] depicting the 
mid-infrared evolution of SN 2004et. The supernova position is 
shown in the pre-explosion SINGS images at day -104 revealing 
evidence of extended emission in this region. The first mid-IR im- 
ages of the SN, obtained 64 days after explosion (second row of 
figure), showed the SN to be very bright. This was during the pho- 
tospheric plateau phase which characterises Type II-P supernovae, 
where hot blackbody emission dominates the optical emission and 
its Rayleigh- Jeans tail extends into the infrared. By day 1015 (third 
row), it can be seen that the SN has faded to almost pre-explosion 
levels, but a late rise in brightness is clearly evident by day 1395. 
Comparable difference images depicting the net mid-IR emission at 
the SN position for the same epochs are also shown in Figure[5] The 
rightmost three boxes in the top row again show the pre-explosion 
images at day -104, which were used as the reference images for 
subtraction from the post-explosion images, to yield the difference 
images shown in the remaining panels. At day 1015 when the SN 
has faded to its faintest levels, a detection at 4.5 pm can just be 
discerned. For the same epoch at 5.8 pm, whilst there is positive 
emission coincident with the position of the SN, this is at a similar 
level to the average noise levels in the residual background of the 
difference image and is therefore considered to be a non-detection 
for which an upper limit to the flux is derived. The SN is much 
more clearly detected in the 8.0-/im difference image at day 1015, 
although the irregular residual background, seen as diagonal bands 



across all of the difference images as this wavelength, creates rela- 
tively large uncertainties in the final measured flux. 

For the 16- pm Spitzer-PUI data, the pre-explosion flux esti- 
mated by Kotak et al. (2009 ) was ^56% higher than that used 
here and so we obtain generally higher post-explosion fluxes than 
those presented by Kota k et ah] p009| ). The factor by which our 
post-explosion fluxes exceed those of |Kotak et al.| ( [2009] > varies 
from ~ 1.1 to 3.3. The largest differences in the photometry occur 
for days 1007-1044, when the SN was close to its faintest levels 
in this waveband. During the final epochs observed (days 1212— 
1173), when the rebrightening was strongest, the fluxes agree to 
within 10%. 



For the 24-pm Spitzer-MlPS data, Kot ak et al.| {2009} mea- 
sured the flux in the two SINGS pre-explosion images observed at 
days -75 and -73, as processed with the standard Spitzer pipeline. 
They found the flux at day -75 to be a factor of almost 1.3 higher 
than that at day -73. From these they measured an average pre- 
explosion flux which is almost 40 % lower than the pre-explosion 
flux measured here. As discussed in Section [AT] the pre-explosion 
MIPS image used here for the analysis of SN 2004et was the 
SINGS enhanced mosaic from their 5th data delivery, which com- 
bined data from both dates. We found a flux consistent with that 
of |Kotak et al. (2009 ) at day -73 but at day -75 we measured the 
flux to be about a factor of 1.7 higher. We also measured the fluxes 
of a number of nearby isolated point sources that were present in 
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Figure 5. Leftmost three columns: the SN position in Spitzer images at selected epochs, showing the mid-IR evolution in IRAC bands 4.5, 5.8 and 8.0/xm. 
The first row shows the SN field in the pre-explosion SINGS images at day -104. The first mid-IR detections were obtained 64 days after explosion when the 
SN was still very bright (second row). The third row shows that the SN had faded to almost pre-explosion levels by day 1015, yet by day 1395 the SN has 
brightened again (fourth row). Rightmost three columns: difference images for the same epochs and wavelengths. The pre-explosion images at day —104 (top 
row) have been registered to, and subtracted from, the post-explosion images using PSF-matched difference imaging techniques (see Section |AT) . Note the 
strong, uneven background residuals at 8.0 fim. 



both images. No systematic offset was found, with the photometry 
differing by no more than 1-10%. 

Despite measuring a higher pre-explosion flux than Kotak 
|et al.| ( |2009| > at 24 /im, we find that photometry from our differ- 
ence images yielded generally higher post-explosion fluxes than 
theirs. Different techniques/software were used to produce the fi- 
nal difference images, which may have resulted in different back- 
ground residuals affecting the photometry. In addition, Kota k et aT] 
(2009 ) used aperture photometry to measure their fluxes, whereas 
we used PSF-fitting. A comparison between results from aperture 
photometry and PSF-fitting is presented by Fabbri (201 1). In addi- 
tion, different techniques may have been used to interpolate the 16- 
and 24-/xm data to the IRAC epochs used for the SED analysis. 



5 OPTICAL AND NEAR-INFRARED PHOTOMETRY 
5.1 The observations 

Optical and NIR photometric observations of SN 2004et were ob- 
tained as part of the SEEDS program over the years 2004 to 2009, 
spanning 79-1803 days after explosion. 

The first two epochs of optical photometry were obtained dur- 
ing the plateau phase of the SN at days 79 and 89 with the 32-inch 
Tenagra II telescope in Arizona. A further three epochs of opti- 
cal photometry were obtained during the nebular phase between 
days 317 and 664 with the Gemini Multi- Object Spectrograph on 



Gemini-North (GMOS-N) in Hawaii. In addition, an archival flux- 
calibrated Subaru spectrum was used to obtain optical photometry 
at day 646 by integrating over the BVRI filter transmission curves. 
Two epochs of late-time, high resolution observations of the SN 
field were obtained with the HST Wide Field Planetary Camera 2 
(WFPC2) around three years after explosion at days 1054 and 1215. 
A final epoch of optical photometry, almost 4 years after explosion, 
was obtained with GMOS-N at day 1412. 

The first epoch of NIR photometry was obtained at day 268 
with the 2.3-m Bok telescope, part of the Steward Observatory at 
Kitt Peak, Arizona. This was followed by three epochs of data from 
the Near InfraRed Imager (NIRI) on Gemini-North at dates close 
in time to the optical images taken with GMOS-N during the neb- 
ular phase. Late-time, high-resolution data was obtained with the 
HST Near Infrared Camera and Multi-Object Spectrometer 2 (NIC- 
MOS2), at epochs corresponding to those of the optical WFPC2 
data. A final i7-band image was taken approximately 5 years (day 
1803) after explosion with the WIYN High-resolution InfraRed 
Camera: WHIRC ( |Meixner et al.|2010) , on the WIYN 3.5-m tele- 
scope at Kitt Peak, Arizona. 

Table [4] provides a complete log of the optical and NIR pho- 
tometric observations of SN 2004et taken as part of the SEEDS 
project. Appendix |A2| describes how the optical and near-IR data 
were processed, while Appendix |A3| describes our late-time HST 
optical and near-IR images which revealed the single point source 
seen at the SN position in IRAC images to be comprised of at least 
three sources. 
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Table 4. Log of optical and near-infrared photometric observations of SN 2004et from the SEEDS collaboration. 



Date 


Age 


Telescope/ 


Filters 


Exp. time 


Program ID 


Principal 




[days] 


instrument 








Investigator 


2004-12-10 


79 


Tenagra II 32" 


VRI 


9x100s 


- 


D. Welch 


2004-12-20 


89 


Tenagra II 32" 


VRI 


9x100s 


- 


D. Welch 


2005-06-17 


268 


Bok 2.3m IR Camera 


JHK 


20x30s 


- 


K. Gordon 


2005-08-05 


317 


Gemini GMOS-N 


g r i 


1x60s 


GN-2005B-Q-54 


G. Clayton 


2005-08-05 


317 


Gemini NIRI 


JHK 


22x30s 


GN-2005B-Q-54 


G. Clayton 


2005-10-17 


390 


Gemini NIRI 


JHK 


10x30s 


GN-2005B-Q-54 


G. Clayton 


2005-10-31 


404 


Gemini GMOS-N 


g r i 


1x60s 


GN-2005B-Q-54 


G. Clayton 


2006-07-06 


652 


Gemini NIRI 


JHK 


22x30s 


GN-2006A-Q-52 


G. Clayton 


2006-07-18 


664 


Gemini GMOS-N 


g r i 


1x60s 


GN-2006A-Q-52 


G. Clayton 


2007-07-08 


1019 


HST WFVC2 


F606W,F814W 


4x400 s 


GOH229 


M. Meixner 


2007-07-08 


1019 


HST N1CMOS2 


F110W,F205W 


5x128s 


GOH229 


M. Meixner 


2007-07-08 


1019 


HST N1CMOS2 


F160W 


4x128s 


GOH229 


M. Meixner 


2008-01-20 


1215 


HST WFPC2 


F606W,F814W 


4x400s 


GOH229 


M. Meixner 


2008-01-20 


1215 


HST NICMOS2 


FU0W,F205W 


5x128s 


GOH229 


M. Meixner 


2008-01-20 


1215 


HST N1CMOS2 


F160W 


4x128s 


GOH229 


M. Meixner 


2008-08-04 


1412 


Gemini GMOS-N 


g r i 


2x600 s 


GN-2008B-Q-44 


M. Barlow 


2009-08-30 


1803 


WI YN WHIRC 


H 


4x180s 


2009B-0516 


M. Otsuka 



5.2 Photometry and light curve evolution 

The final optical VRI magnitudes for SN 2004et are presented 
in Table [5] including those of neighbouring star 2 discussed previ- 
ously. The magnitudes for epochs from day 646 onwards have been 
corrected for the contribution from star 2, but this was not necessary 
for earlier epochs when the SN brightness dominated. The B-band 
magnitude at day 646, corrected for the contribution from star 2, is 
detailed in the notes to the table. The NIR JHK magnitudes of the 
supernova are given in Table [6] although the magnitudes measured 
from the NICMOS images at days 1019 and 1215 are HST Vegam- 
ags (Appendix A2). The optical and NIR light curves are presented 
in Figures [6[and[7]respectively. 

Figure|6] combines the BVRI data of |Sahu et al] ( |2006| » with 
the SEEDS data described above to provide optical light curves 
(open symbols) from a few days after the explosion to almost 4 
years later. The light curves are well sampled until around day 400. 
Few Type II SNe have been observed beyond ~ 500 days. However, 
SN 1987 A, the closest supernova to have occurred in the past cen- 
tury, has been well studied for over 2 decades and provides detailed 
light curves f or comparison. The broad-band BVRI light curves 
of SN 1987A ( |Hamuy & Suntzeff|l990l[Walker & Suntzeff|l99l| 2 | 
are plotted as solid curves in Figure [6] normalised to the data of 
SN 2004et at around 200 days. The light curve evolution of both 
supernovae is quite similar from the early nebular phase (~ 160 
days) to around day 650, although beyond this time their evolution 
is markedly different. The brightness of SN 2004et clearly levels 
off from around day 1000 in the VRI bands (the last B-band mea- 
surement was at day 646), while SN 19 87 A continued to fade. In- 
terestingly, the slope of the radioactive decay curves resemble those 
of SN 1987 A between days ~ 720 and 800, around the time when 
dust production for SN 1987 A was assumed to have ended (day 
775; |Wooden et al.|1993) . 

For the first few years during the nebular phase, the light curve 



of Type II SNe is predominantly powered by 7-rays from the ra- 
dioactive decay of 56 Co to 56 Fe, at a rate corresponding to the e- 
folding time of the 56 Co decay (T56 = 111.3 days). For example, 
the i?-band photometry of the Type II SN 1990E ( |Benetti et al] 
1994 ) closely follows this evolution through to ^540 days post- 
explosion, suggesting that simple 56 Co decay provides a good esti- 
mate of the unextinguished i?-band light curve for at least that long. 
The expected decay rate is 7 (mag per 100 days) = 0.98 for com- 
plete 7-ray trapping ( |Patat et al.|1994| >. For SN 2004et, [Sahu et al.| 
('2006) found that the decline of the broad-band BVRI magnitudes 
during the early nebular phase (180-310 days) was linear, with 
decay rates of 7^ ~0.64, ~ 1-04, 7^ ~ 1.01 and 7/ ~ 1.07. 
Maguire et al.| p010| > found similar results from their own data 
(~ 136-300 days), with lB = 0.64 ± 0.02, j v = 1.02 ± 0.01, 
7h = 0.92 ± 0.01 and 71 = 1.09 ± 0.01. With the exception 
of the B band, the decay rates were close to that of 56 Co decay, 
suggesting that 7-ray trapping was efficient during this phase. 

However, both authors found that the optical decay rates of 
SN 2004et steepened beyond ~ 300 days, suggesting that either the 
supernova had become transparent to 7-rays and hence 7-rays as a 
source of energy were escaping, or that dust was forming within 
the ejecta and causing localised optical extinction, or possibly was 
due to both phenomena. 

To investigate this further we can look again at the radioactive 
deposition. As the ejecta expand, their opacity to 7-rays is expected 
to decrease, which results in a modified light curve of the form 
fWoosley etal.|1989| >: 



LIS) 



^56 4>0 (to/t) 



(1) 



2 |Hamuy & Suntzeff| l 1990) data downloaded from the NOAO archive: 
ftp://ftp.noao.edu/snl987a/ubvri.txt 



where the term in the brackets is the deposition function, i.e., the 
fraction of 7-rays deposited in the envelope; ^56,7 = 0.033 cm 2 g _1 
is the average opacity to 56 Co-decay 7-rays, and 0o = 7xl0 4 g 
cm -2 is the column depth at the fiducial time to = 1 1.6 days chosen 
to match the bolometric light curve of SN 19 87 A. 

After this first source of decay energy has become sufficiently 
weak, other energy sources which could become important in pow- 
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Table 5. Optical photometry of SN 2004et. 



UT date 


Age 




Magnitudes 




Source 




[days] 


V 


R c 


Ic 




2004-12-10 


79 


13.09 ± 0.04 


12.38 db 0.03 


11.93 ±0.03 


Tenagra II 


2004-12-20 


89 


13.21 ± 0.03 


12.45 ± 0.03 


11.98 ±0.03 


Tenagra II 


2005-08-05 


317 


17.35 ± 0.02 


16.52 db 0.03 


15.91 ± 0.04 


Gemini GMOS-N 


2005-10-31 


404 


18.28 ±0.04 


17.69 ± 0.03 


16.87 ± 0.04 


Gemini GMOS-N 


2006-06-30 


646* 


21.59 ±0.55t 


21.00 ±0.551" 


20.63 ± 0.56t 


Subaru FOCAS 


2006-07-18 


664 


22.13 ±0.06t 


21.56 ±0.09* 


21.37 ±0.1lt 


Gemini GMOS-N 


2007-07-08 


1019 


23.20 ± 0.20* 




22.70 ± 0.20* 


HST WFVC2 


2008-01-20 


1215 


23.40 ± 0.30* 




23.00 ± 0.30* 


HST WFPC2 


2008-08-04 


1412 


23.80 ± 0.40t 


22.87 ± 0.281" 


22.80 ± 0.70* 


Gemini GMOS-N 


Star 2 




24.2 ± 0.3 


23.5 ± 0.5 


22.9 ± 0.4 


HST WFVC2 



Evans < 1993 »; Glass 1 1999 



Zero-magnitude flux [Jy] 3670.3 2972.3 2402.1 

X ef f [H 055 064 O80 

* Optical photometry at day 646 was estimated from an archival Subaru-FOCAS spectrum by integrating over the BV RI filter transmis- 
sion curves. The B band magnitude of SN 2004et at this time was 22.47 ± 0.22, corrected for an estimated S-band contribution from 
star 2 of 25.46 ± 0.50 (see text). 

1" The SN magnitudes at these late epochs have been corrected for contamination by star 2 (whose magnitudes measured from the high- 
resolution HST data are listed in the final table entry). 

$ Since the SN and star 2 were resolved in the day 1019 and day 1215 HST images, the magnitudes given for these epochs are for the SN 
alone. 
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Figure 6. BVRI light curves of SN 2004et, based on the data of Sahu et al. (2006) and the SEEDS photometry listed in Table[5] For clarity, the plots have 
been vertically shifted by the amounts shown in the legend. The black dashed lines are the expected light curves based on the decay of 56 Co and other isotopes 
jWoosley et al.|1989) normalised to the V and R band magnitudes of SN 2004et during the early nebular phase (between ~ 160-260 days). The solid lines 
are the corresponding light curves of SN 1987A jHamuy & Suntzeff|l 990] | Walker & Suntzeff|199l) normalised to those of SN 2004et at ~ 200 days. Dotted 
lines are only to guide the eye along the horizontal or "flat" parts of the light curve. These curves are typical for every SN that shows a light echo, in that the 
light curves level off or "flatten" once the SN flux drops below that of the echo. 



ering the very late time light curves are 7-rays, positrons and elec- 
trons from the radioactive decay of 57 Co, 44 Ti and 22 Na. The equa- 
tions that describe the energie s from all thes e isotopes, including 

following the work 



1993 



Co, were summarised by |Li et al. 
of Woosley et al. (1989 ), to describe the deposition behaviour of 
SN 19 87 A. Adopting the same deposition behaviour for SN 2004et, 
the radioactive decay curve attributable to the energy sources from 
these isotopes, including a term to account for the decrease in opac- 
ity to 7-rays as the ejecta expands, has been plotted in Figure [6] 



(dashed line) over both the V and R band magnitudes of SN 2004et, 
normalising to the early nebular phase data (^160-260 days). As 
expected, with the exception of the B band, the decay rates dur- 
ing the early nebular phase closely follow those of the radioac- 
tive decay deposition. However, from about 400 days the R band 
light curve has clearly begun to decline more rapidly than the ex- 
pected light curve from radioactive decay deposition. The / band 
follows a similar trend, whereas the steepening of the decline rate 
appears to occur slightly later in the V band, having clearly be- 
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Table 6. Near-infrared photometry of SN 2004et. 



UT date 


Age 




Magnitudes 




Source 




[days] 


J 


H 


K 




2005-06-17 


268 


15.14 ±0.05 


15.15 ±0.03 


... a 


Steward/Bok IR Camera 


2005-08-05 


317 


16.01 ± 0.03 


15.84 ± 0.03 


15.18 ±0.04 


Gemini NIRI 


2005-10-17 


390 


16.96 ± 0.03 


16.62 ± 0.03 


16.23 ± 0.04 


Gemini NIRI 


2006-07-06 


652 


20.09 ± 0.05 


19.52 ± 0.06 


19.19 ±0.07 


Gemini NIRI 


2007-07-08 


1019 


22.25 ±0.13 b 


22.61 ± 0.36 b 


21.91 ±0.28 b 


HST NICMOS 


2008-01-20 


1215 


22.55 ±0.14 b 


22.69 ± 0.43 b 


21.44 ±0.18 b 


HST NICMOS 


2009-08-30 


1803 




<22.6 




WIYN WHIRC 


Zero-magnitude flux [Jy] 


1656.3 


1070.9 


672.8 


[Glass] (l999| 


X eff [/xm] 




1.25 


1.65 


2.20 



a There is no K band magnitude at day 268 as unusual image artifacts compromised the photometry. 

b HST Vegamags in NICMOS 2 filters F110W (~ J band), F160W (~ H band) and F205W (~ K band) for days 1019 
and 1215. For each filter and epoch, the measured count rate (CR, in units of DNs - 1 ) at the position of the SN was converted to 
flux by multiplication with the PHOTFNU (Jy sDN — 1 ) conversion factor given in the fits header, where PHOTFNU is 
the bandpass-averaged flux density for a source that would produce a count rate of 1DN — 1 . PHOTFNU = 1.21 x 10 — 6 , 
1.50 X 10 -6 and 9. 69 X 10 -7 JysDN -1 for F110W, F160W and F205W respectively. 
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Figure 7. Late-time NIR light curves of SN 2004et, based on the SEEDS photometry listed in Table[6] For clarity, the plots have been vertically shifted by the 
amounts shown in the legend. The solid lines are the corresponding light curves of SN 1987A (Bouchet & Danziger 1993 1 normalised to the light curves of 
SN 2004et at day 317, during the early nebular phase. Dotted lines are only to guide the eye along the horizontal or "flat" parts of the light curve. The inset 
shows the if -band light curve only for the extended period to day 1803, the final SEEDS NIR observation, taken with the WIYN telescope in August 2009. 
Contamination from neighbouring sources observed in the high-resolution NICMOS images at days 1019 and 1215 resulted in an ambiguous detection of the 
SN with the WHRIC at day 1803, for which an upper limit (downward pointing arrow) has been derived. 



gun sometime between 460 and 540 days. In comparison with the 
expected radioactive decay deposition behaviour of SN 1987 A, as 
modelled by |Li et aE| ( |1993| ), there is evidence for a steepening de- 
cline of the light curves, indicative of dust formation in the ejecta of 
SN 2004et from around 400 days, and possibly earlier (Sahu et al. 
|2006[|Maguire et al.|20i0| . From the F-band light curve (Figure[6}, 
we estimate that the difference between the observed and predicted 
light curves was 0.8 magnitudes by day 690, if allowance is made 
for the effective opacity term for 56 Co 7-rays, or 1.5 magnitudes if 
this term is neglected. 

It is clear that the elevated brightness of SN 2004et after 1000 
days cannot be explained by the inclusion of isotope decays, such 
as 57 Co and 44 Ti, which could be important at these late times. The 



plateauing of the optical light curves above the expected radioactive 
decay suggests an additional energy source has come into play by 
at least day 1000 and is consistent with the late rise observed in the 
mid-IR observations after this time. Similar phenomena observed 
for other SNe have been attributed to light echoes (e.g., SN 2007od; 
|Andrews et al.|2010| >. 

Figure [7] presents the late-time NIR light curves (open sym- 
bols) of SN 2004et based on the SEEDS data described previ- 
ously. The JHK light curves of SN 1987A (solid lines; |Bouchet| 
|& Da nziger 1993]) have been arbitrarily scaled to the early nebu- 
lar phase data of SN 2004et for comparison. The NIR light curves 
of SN 2004et are not well-sampled but clearly deviate from those 
of SN 1987 A by day 646, after which time SN 2004et is system- 
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atically brighter. The NIR light curve evolution reflects that of the 
optical, with a relative plateau in brightness occurring from around 
1000 days, consistent with a light echo hypothesis. 

However, by day 1803 (see inset in Figure [7]), the if -band 
brightness has faded to beyond a clear detection with the WHRIC 
detector on the 3.5-m WIYN telescope. The derived upper limit of 
22.6 mag accounts for contamination by neighbouring stars (as de- 
scribed in Section |A2| ) and suggests the SN faded sometime after 
day 1215. 



6 SPECTRAL ENERGY DISTRIBUTION ANALYSIS 
BASED ON OPTICAL, NIR AND MID-IR 
PHOTOMETRY 

The mid-infrared photometry listed in Table [3] along with the opti- 
cal and IR photometry in Table [5] and Tableland the optical pho- 
tometry of Sahu et al. ( 2006 ), have been used to investigate the 
spectral energy distribution of SN 2004et at each of the IRAC ob- 
servation epochs from days 64 to 1395. 



6.1 Blackbody fitting 

To investigate the physical processes that determine the observed 
optical and infrared continuum emission and their evolution, black- 
bodies were matched to the SEDs at each of the IRAC epochs from 
day 64 to day 1395. Where necessary, the light curves for the op- 
tical, NIR, PUI 16-/xm and MIPS 24-/mi data were used to in- 
terpolate their measured flux densities to the epochs of the IRAC 
data. It should be noted that to extrapolate the last B magnitude, 
obtained on day 646 to the closest IRAC epoch at day 690, the 
better- sampled F-band decline rate during this period was adopted 
in order to account for the gradual flattening of the light curve. Sim- 
ilarly, the JHK magnitudes from day 652 were extrapolated to the 
closest IRAC epoch of day 690 assuming the decline rate observed 
in the 7 c -band during this period. Gemini Michelle flux densities at 
11.2 fim were not interpolated due to insufficient data but, where 
available, they are compared with the closest IRAC epoch. 

The interpolated optical data were converted from the stan- 
dard Johnson-Cousins BVR C I C magnitudes to flux densities us- 
ing the zero-magnitude flux densities of Evans (1993) and Glass 
(1999). The interpolated Steward and Gemini NIR data were 
converted from standard JHK magnitudes to flux densities us- 
ing the zero-magnitude flux densities of |Glass| ([T999). The late- 
time ffiT-NICMOS flux densities were obtained by multiply- 
ing the count rate measured for the SN by the PHOTFLAM 
(erg cm -2 A -1 DN -1 ) conversion factor from the fits image head- 
ers, where PHOTFLAM is the bandpass-averaged flux density 
in F\ for a source that would produce a count rate of 1 DN s _1 . All 



flux densities were de-reddened using E(B-V) = 0.41 mag ( Zwit- 
ter et al. 2004) and assuming the extinction law of Cardelli et al. 
( |1989| > with R v = 3.1, corresponding to A v = 1.27 ± 0.22 mag. 



6.1.1 Day 64 

Panel (a) of Figure [8] shows a fit to the SED at day 64, during 
the photospheric plateau phase of SN 2004et. The fit, which uses 
a 5400 K, blackbody normalised to the IRAC 5.8-/xm flux den- 
sity, under-estimates the V and R band flux densities, which we 
attribute to the presence of strong e mission lines such as Ha in 
the 5400-7000 A wavelength region ( |Sahu et al.|2006||Kotak et al. 



2009), while it significantly overestimates the U band flux den- 
sity. A spline curve fitted to the UB photometry was therefore used 
to estimate the total flux shortwards of B. This spline curve was 
combined with the 5400 K blackbody, truncated at wavelengths 
< B (4400 A), to yield a total integrated flux of 4.5 x 10" 13 W m" 2 
(with a corresponding luminosity of 4.9xlO 8 L ), about 93% of 
the total integrated flux/luminosity from the hot blackbody alone. 
At wavelengths < B, the total flux corresponding to the spline fit 
constituted about 49% of that from the 5400 K blackbody fit. In 
panel (a) of Figure [8] the solid line represents the combined spline 
and truncated blackbody fit, whilst Table [7] lists the parameters for 
the blackbody fits. 

It is well known that all Type I SNe show a pronounced early- 
time deficit at ultraviolet (UV) wavelengths relative to a blackbody 
fitted at longer wavelengths (e.g., [Panagia 2003). This has been in- 
terpreted as due to strong line blanketing by many low excitation 
lines of Fell and other lines shortwards of ~ 4000 A (e.g., Branch 
|& Ve nkatakrishna 1986). The situation for Type II SNe seems to be 
less clear. In his review of optical spectra of supernovae, Filippenko 
( |1997| ) noted that most Type II SNe do not show this feature, with 
the early-time spectra approximating a single-temperature Planck 
function from UV through to IR wavelengths, and occasionally 
even showing a slight UV excess. However, Franss oiTet al.| (l987), 
from their studies of the peculiar Type II SN 1987A, concluded 
that supernova atmospheres with a normal (solar) chemical compo- 
sition can give rise to line blanketing effects, such as those seen in 
the UV spectra of SN 19 87 A. They proposed that the differences 
in UV spectra of supernovae may instead be due to differences in 
the density of the CSM. The earliest IUE (1 150-3200 A) spectra of 
SN 1987 A (rsj 6 days after explosion) showed a strong UV deficit in 
the wavelength range ~ 1250-3200 A, relative to the 6000 K black- 
body curve defined at optical and infrared wavelengths ( Danziger 
|et al.|1987] >. This was still present at day 60, and possibly as late 
as day 260, as shown by the best-fit SEDs of |Wooden eTaL] ( [T993] ). 
The apparent UV drop-off relative to the blackbody fit to the photo- 
metric SED of SN 2004et at day 64 suggests that a similar effect is 
present in this Type II-P SN. |Li et aL| ( |2005| ) noted from a spectrum 
of SN 2004et at day 9 that "there is a peculiar decline blueward of 
4000 A not commonly observed in the spectra of normal SNe II-P". 



6.1.2 Days 300-828 

The blackbody fits to the SEDs at days 300 to 828 are shown in 
Figure[8] The corresponding blackbody parameters are listed in Ta- 
ble [7] For each epoch, more than one component was required to 
obtain a reasonable fit to the optical, NIR and mid-IR photometry. 
These were comprised of (i) a hot blackbody, with temperatures 
during the period 300-828 days ranging from 7100-10000 K, rep- 
resenting the optical and NIR continuum emission from the op- 
tically thick hot gas of the ejecta, (ii) a warm blackbody (420- 
650 K) representing the emission at mid-infrared wavelengths, and 
(iii) from day 464 onwards, a cooler (200-250 K) blackbody to ac- 
count for the observed emission at the two longest mid-IR wave- 
lengths (16 and 24 /mi). Excess emission is also present at 24 /am 
on days 300-406, but at a relatively constant level, see Figure|8] 

From days 300 to 464, the mid-IR emission demonstrated a 
clear excess at 4.5 /xm, as noted by Kotak et al. (2009). A simi- 
lar feature has been seen in mid-IR photometry of other Type II 
SNe, such as the Type II-P SN 2 003gd at day 499 (|Sugerman et al.| 
12006] ), and SN 2007it at day 340 ( [Andrews et al.|2010) . This can be 
attributed to emission from the carbon monoxide (CO) fundamen- 
tal band at 4.65 /mi, which was directly observed in the infrared 
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Wavelength [/xm] Wavelength [/xm] 

Figure 8. Blackbody fits to the day 64 to 1395 SEDs of SN 2004et as defined by the optical-IR observations (wavelength range plotted: 0.3-30 /im). Filled 
circles indicate fluxes observed at the epoch of the IRAC observations. Open circles indicate optical, NIR and mid-IR photometry which has been interpolated 
or extrapolated to the epochs of the IRAC observations. Filled squares indicate the closest epoch {i.e., not interpolated) HST-NICMOS photometry and filled 
triangles indicate the closest epoch Gemini-Michelle N'-band photometry (see Table [2] for epochs of the Michelle and NICMOS observations). Upper flux 
limits for non-detections are indicated by the downward-pointing arrows. Where error bars are not shown, uncertainties are smaller than the symbol size. All 
flux densities were de-reddened using E(B- V) = 0.41 mag ( Zwitter et al. 2004) and the extinction law of Cardelli et al. (1989) with Ry = 3.1. 



spectra of the Type Ilpec SN 1987 A from as early as 100 days af- 
ter explosion (Suntzeff & Bouchet 1990) and stayed visible until 
at least day 615 (Wooden et al.|1993| ). The red wing of the same 
emission line was detected in the Spitzer-IRS spectra of the Type 
II-P SN 2004dj at days 109 and 129 ( |Kotak et al.|2005] >. Together 
with the detection of the first overtone of CO at ~ 2.3 (im in the 
NIR spectrum of SN 2004et ( [Maguire et al.|2010|), as wel l as in the 
spectra of several other Type II SNe ( Gerardy et al. 2002 and refer- 
ences therein), these observations indicate that strong CO emission 
is common in Type II SNe. 

Kotak et al. ( 2009) also noted an excess at 8.0 /xm compared to 



blackbodies matched to the mid-IR continuum emission from days 
300 to 690. They found this to be consistent with a broad emission 
feature between 8-14 /xm seen clearly in Spitzer-IRS spectra until 
at least days 450-481 and attributed this to silicate emission, with a 
contribution from the silicon oxide (SiO) fundamental band in the 
7.7-9. 5-/xm region. Accordingly, our warm blackbody component 
was normalised to the 3.6 or 5.8-/xm fluxes as best representing 
the continuum emission during these times. For the best-fit SED 
at day 300 (panel (b) in Figure [8}, the Gemini 11.2-//m flux from 
day 3 1 1 is also underestimated by the warm blackbody component. 
This excess is attributed to the broad silicate emission feature seen 
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Figure 9. The temperature evolution of the blackbody components that 
were fit to the SED of SN 2004et at different epochs. Hot component tem- 
peratures ranged from 5400 K to 10000 K, warm component temperatures 
ranged from 370 K to 650 K and cool component temperatures ranged from 
125 K to 250 K (see Table 7). 



at its strongest in the IRS spectra at days 294 and 349 presented by 
Kotak et al. (2009). The 8.0-/xm excess is no longer evident after 
day 690, consistent with the IRS spectrum at day 823 presented by 
|Kotaket al.| ( [2009l >. Panel (g) of Figure[8]shows that by day 828 the 
observed 8.0-/xm flux is well-matched by the blackbody fits. 



6.1.3 Days 1015-1395 

Best fits to the SEDs at epochs of days 1015, 1054, 1 191 and 1395 
are also shown in Figure [8] and the corresponding blackbody pa- 
rameters are listed in Table [7] Hot (10000 K), warm (350-520 K) 
and cool (120-170 K) blackbody components were matched to the 
optical, NIR and mid-IR photometry. For days 1015 and 1191, the 
adopted NIR fluxes were those measured in the F110W, F160W and 
F205W NICMOS2 filters at the reasonably contemporary epochs of 
days 1019 and 1215. 



6.2 Discussion of results from blackbody fitting 

The complete set of parameters from the blackbody fits to the pho- 
tometry from days 64 to 1395 are listed in Table[7] The temperature 
evolution of the best fit multi-component blackbodies is shown in 
Figure [9] 

The left-hand panel of Figure [10| shows the evolution of the 
luminosities of each of the blackbody components, and their sum, 
compared with the theoretical luminosity due to the radioactive de- 
position of 56 Co, 57 Co, and other isotopes, adopting the deposi- 
tion behaviour of SN 1987A as modelled by |Li et al.|(T993] > and 
Woosley et al. (1989). The radioactive decay deposition curve for 
SN 19 87 A was scaled by a factor of 0.69 to normalise it to the total 
luminosity of SN 2004et between days 300-464. It should be noted 
that the blackbody fits to the photometry mainly trace the contin- 
uum emission of the SN (known optical and mid-IR emission line 
features were deliberately not matched by the blackbodies as dis- 
cussed in the previous sections). As a consequence, the luminosities 
listed in Table [7] slightly underestimate the total luminosities. 



6.2.1 Evolution of the hot component 

The day 64 photometry can be fitted by a 5400 K blackbody, con- 
sistent with an origin from the optically-thick emission of the ejecta 
photosphere, with the Rayleigh- Jeans tail of the corresponding 
blackbody extending into the infrared. As the ejected envelope adi- 
abatically expands and cools, the hydrogen ionised by the initial SN 
shock approaches the temperature for recombination (~ 5000 K) 
and a recombination wave recedes through the envelope. The re- 
combination front defines the photosphere and as such the temper- 
ature of the photosphere is characterised by the recombination tem- 
perature of hydrogen. Wooden et al.| ( [T993] ) found that a hot com- 
ponent of 5000-5500 K was a good match to the SED of SN 1987 A 
from days 60 to 777, which they found to be in agreement with ob- 
servations of other Type II supernovae by Kirsh neTet al.|(T973] >. 
With a similar temperature, the 5400 K blackbody fit shown in 
panel (a) of Figure [8] is representative of the photospheric con- 
tinuum emission at this time. The ejecta velocity of 3167kms _1 
implied by the emitting radius of the 5400 K blackbody is in rea- 
sonable agreement with the velocity of ~ 3500 km s _1 found by 
|Sahu et al.] (2006) for SN 2004et in the plateau phase, estimated 
from the minimum of weak, unblended absorption lines of Fe II at 
4924, 5018 and 5169 A. [Kotak et al.||2009l obtained a reasonable 
fit to the day 64 SED of SN 2004et with a single 5300 K blackbody, 
concluding there was little sign of thermal emission from dust. 

As expected, the corresponding luminosity at day 64 exceeds 
that from radioactive decay deposition, which only begins to dom- 
inate the light curve during the nebular phase (from ~ 130 days, 
Maguire et al. 2010), following the sharp decline from the plateau 
at ~ 110 days (Sahu et al.|2006| >. For the 5400 K blackbody fit the 
day 64 luminosity exceeds that from radioactive deposition by a 
factor of 4.9, or by a factor of 4.5 if considering the spline plus 
truncated blackbody fit which is a better match to the U band data. 
This compares to the factor of 3.8 found by Kotak et al. (2009). 

By day 300, the estimated temperature of the hot blackbody 
component had increased to 7250 K and then remained relatively 
constant (± ~ 150 K) until day 464 (Figure|9|, whilst its luminosity 
faded quite rapidly (by a factor of ~ 7) during this time (FigureflO). 
Note that at day 300, the blackbody fit significantly over-estimates 
the fluxes at the shortest wavelengths in the U and B bands, as also 
seen at day 64 (see previous discussion). By day 690 the temper- 
ature had increased again to 8500 K and to 10,000 K by day 828, 
whilst the luminosity continued to decrease (by a factor of ~ 15). 
For epochs beyond 1000 days, the temperature of the hot compo- 
nent appears to remain constant. 



6.2.2 Evolution of the warm and cool dust components 

The warm dust component cooled monotonically from 650-420 K 
between days 300 and 690 (Figure [9]). This is consistent with the 
fading of SN 2004et observed in the mid-IR during this time. The 
sum of the luminosities of the hot component and warm compo- 
nent is less than or comparable to the predicted radioactive depo- 
sition up to day 690. On day 828 however, the sum of the hot and 
warm component luminosities was a factor of 1 .7 higher than the 
predicted radioactive deposition luminosity. 

The presence of a mid-IR excess from day 300, demonstrated 
by the requirement of a warm component to match the SEDs from 
this time, and the evolution of this warm component from 300 to 
690 days, are consistent with emission from dust freshly synthe- 
sised in the supernova ejecta, in agreement with the results of Kotak 
|et al.| ([2009) from their interpretation of similar data. 
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Table 7. Blackbody-fitted parameters to the SEDs of SN 2004et. 1 
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Up to three components - hot, warm and cool - were used to fit the continuum emission of the SEDs (see text). For each component, R is the blackbody radius corresponding to the best fit temperature, T; v is 

the expansion velocity corresponding to the radius R; and L is the luminosity for an adopted distance to the SN of 5.9 Mpc. 
1 The day 64 5400 K blackbody fit over-estimated the U band flux. A spline curve fit to the UB fluxes and extrapolated to shorter wavelengths was combined with the 5400 K blackbody fit to the longer wavelength 

data, where the blackbody was truncated at wavelengths shorter than the £?-band - see panel (a) of Figure[8] The total integrated flux and corresponding luminosity of this "spline + 5400 K blackbody" fit are 

about 93% of those values for the pure 5400 K blackbody fit shown in this table. 
3 At day 300, it was found that the hot component blackbody still over-estimated the fluxes at the shortest wavelengths. A spline curve fit to the UBV fluxes and extrapolated to shorter wavelengths was combined 

with the two-component blackbody fit to the longer wavelength data, as shown in panel (b) of Figure[8] where the 7250K blackbody component was truncated at wavelengths shorter than the V-band. The total 

integrated flux and corresponding luminosity of this "spline + 7250K + 650K blackbody" fit are about 73% of those values for the combined 7250K + 650K blackbody fit shown in this table. 
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Figure 10. The luminosity evolution of SN 2004et, from multi-component blackbody fits to the SED at each epoch (see Table 7). The temperature evolution 
of the individual components is shown in Fig. 9. Left: The luminosity evolution of the individual components (hot, warm and cool) is shown at each epoch, 
together with that of the corresponding total luminosity. The solid curve shows the predicted luminosity due to radioactive deposition, as scaled from the 
models of |Li et al7|(1993) (see text for details). Right: the sum at each epoch of the luminosities of the two components, warm and cool, attributed to dust 
emission, showing a steady decline up to day 828, followed by a subsequent rise. 



The approximately constant luminosity of the warm compo- 
nent between days 828 and 1054 corresponded to the lowest level 
reached, since at day 1191 the warm component luminosity had in- 
creased by a factor of 4.3. From day 1015 the warm component 
luminosity clearly exceeds that of the radioactive deposition (by a 
factor of ~ 8) and at day 1191 by a factor of over 100. 

The cool dust component, required to fit the mid-IR fluxes 
long wards of 16 /im from day 464 onwards, showed a monotonic 
decline in temperature from 250 K at day 464 to ~ 120 K at day 
1015, but then increased slightly and remained around 160K be- 
tween days 1054 and 1395. Its luminosity stayed roughly constant 
from days 464 to 1054, but had increased by a factor of ~ 2 by days 
1191-1395. The velocities implied by the minimum emitting radii 
of these cooler blackbodies ranged from ~ 2500-6000 km s _1 . 

Kota k et al.] (2009), on the other hand, found that the tem- 



perature of their cool component remained approximately con- 
stant, at 120±10K for all epochs from day 300-1395, with 
minimum blackbody radii corresponding to velocities as large as 
12000 kms -1 , so they ruled out ejecta-condensed dust as a source 
of the cool dust emission. The differences between our results for 
the cool component and those of Kotak et al. ( 2009) are most likely 
due to the differences in the photometry at the longest wavelengths, 
16 and 24 /xm, as discussed in Section |4~2| Our 16 and 24-/xm fluxes 
required generally higher temperature blackbodies to fit them than 
did those of |Kotak e t al. (2009). Since our luminosities for the 
cool component are generally consistent with those of Kota k et aL] 
( 2009), the higher blackbody temperatures that we obtain resulted 
in lower minimum radii and therefore lower minimum outflow ve- 
locities. 

The right-hand panel of Figure [l0| shows the evolution of the 
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Table 8. Input parameters for R 2 density distribution dust models. 
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8000 
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3.5 


1.30 
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sum of the luminosities of the two dust components (warm plus 
cool). It shows a steady decline until day 828, a period of constant 
luminosity until day 1054 and a steep rise thereafter. The evolu- 
tion beyond day 828 implies the presence of an additional energy 
source. Kotak et al. (2009) found similar results and interpreted the 
late rise in the mid-IR flux as due to ejecta-CSM interaction and 
the subsequent formation of a cool dense shell of dust behind the 
reverse shock of the supernova. An alternative explanation in terms 
of a light echo from pre-existing extended CSM dust will be pre- 
sented by Sugerman et al. (in preparation). 



7 RADIATIVE TRANSFER MODELLING: DAYS 300-690 



Table 9. The contribution from line emission to the Vi?/-bands. 



In agreement with previous investigators ( |Sahu et al. 2006 ; Kotak 
|et al.||2009| >, the development of a red-blue asymmetry in optical 
emission line profiles, the evidence from optical light curves for 
additional extinction by dust in the ejecta, and the development of 
a mid-IR excess attributable to dust emission all support the infer- 
ence that dust formed within the ejecta of SN 2004et from about 
300 days after explosion. As discussed above, by day 828 the to- 
tal luminosity exceeded the estimated radioactive deposition lumi- 
nosity by a factor of four (Figure [ToJ, indicating that an additional 
component dominated by that date. Consequently, epochs later than 
day 690, when the thermal IR emission can no longer be solely at- 
tributed to internally heated newly-formed dust in the ejecta, will 
not be modeled here - Sugerman et al. (in preparation) will present 
a circumstellar light echo model for these later epochs. 

To investigate the time evolution of dust formation in the 
ejecta of SN 2004et and to estimate the mass of dust present, a 
number of dust shell models were built to match the observed SEDs 
at the epochs between days 300 and 690 (Table [8}. The models 
were calculated using the three-dimensional Monte Carlo radiative 
transfer (RT) code MOCASSIN prcolano et al.|2003||2005) which 
accounts for the primary and secondary components of the radia- 
tion field in a fully self-consistent manner through absorption, re- 
emission and scattering of photons. The photon paths are followed 
from a specified source through a given composition, grain-size dis- 
tribution, density and geometry of dust. The particular choices of 
these parameters are either constrained a priori or are varied until 
the model emission and extinction match the observations. 

For the day 300-690 models we assumed that the observed 
IR emission originated from dust formed in the SN ejecta. Heat- 
ing is due to 7-rays from the decay of 56 Co, which are reprocessed 
to optical and UV wavelengths through interaction with the gas. 
We assume that this leads to a local radiation field whose strength 
is proportional to the local ejecta density. Based on these assump- 
tions, and following the previous modelling of SN 2003gd by |Sug-| 
|erman et al.| ( [2006| ) and of SN 1987A by |Ercolano eTaL] {2007^ 7 
the RT models were constructed such that the dust and source lu- 
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minosity were mixed within a spherical expanding shell of inner 
radius Ri n and outer radius R ou t = YR in . For the dust density 
distribution, we considered two cases: (i) a smooth distribution fol- 
lowing an r~ 2 density profile (smooth model), and (ii) a clumpy 
model where dense clumps exist in a less-dense interclump medium 
(ICM), where the ICM follows a smooth r~ 2 density distribution, 
with the local heating source located only in the ICM. For each 
case, we compared the observed SEDs with those reproduced by 
the models to determine which gave the best fit to the observations. 

Before running the models, we estimated the contribution 
from line emission to the VRI -bands. We used day 417 and 454 
optical spectra downloaded from SUSPECT plus day 336 TNG 
and day 646 Subaru spectra, along with the Subaru/FOCAS V Ri- 
band filter transmission curved because these have similar band 
centres and widths to the standard Johnson filters. The percentage 
line contributions in each band are listed in Table [9] When evalu- 
ating the fitting accuracy of the SED modeling, we excluded the R 
and / bands on days 300, 360, 406, and 464 because both bands 
were dominated by line emission at those epochs. We also omitted 
the IRAC 4.5-/xm data-points from the fitting, due to the potentially 
very large CO line emission contributions in that band ( Kota k et aT] 
2009 ). Although our modelling took into account potential emis- 
sion from the broad silicate 10-/xm band, the SiO fundamental vi- 
brational band can also contribute to the IRAC 8-/xm band - from 
Spitzer IRS spectra its contribution was deduced by |Kotak et"aL] 
(2009) to be significant on days 300-464. 



7.1 Smooth dust distribution models 

Both amorphous carbon (AC) and silicate dust grains were con- 
sidered. Optical constants were taken from |Zubko et al.|(|1996|) for 
amorphous carbon (their ACH2), and from Drain e & Lee| {|1984 ) 
for the silicates. To investigate the dust composition, models were 
run with amorphous carbon: silicate mixtures of (100-x):x%, for 
x=0,20,40,60,80,100. It was found that the case of 20 % amorphous 
carbon and 80 % silicate (by mass) best matched the observed SEDs 
at all epochs, and this composition was adopted for all subsequent 

models. _ 

A standard MRN a -3 5 distribution ( [Mathis et al.|l977| l with 



See http://www.naoj.org/Observing/Instruments/FOCAS/camera/filters.htrnl 
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a m m=0.005 fim and a max =0.25 /xm could not reproduce the steep- 
ness of the observed SED in the JH if -bands before day 690, nor 
the 16- and 24-/xm flux densities. An improved fit was found using 
an MRN distribution with a m i n =0.1 /im and a max =1.0 /xm and this 
was subsequently adopted for all epochs. 

The density distribution in SN ejecta can range from approx- 
imately flat, to very steep in the layers that had formed the photo- 
sphere of the progenitor. Here, we adopt an r~ 2 law. We found that 
distributions steeper than r~ 3 led to too much emission in the 1.6- 
3.6-/xm spectral region at the earlier epochs, as a result of the high 
densities and high heating rates at the inner edge of the ejecta. As 
a first guess, we adopted initial values for Ri n from our blackbody 
fitting, and then varied them to match the observations. We calcu- 
lated models with shell size scaling parameters Y = R ou t / Rin = 
3.0, 3.2, 3.5, 4.0, 4.5. The adopted heating luminosities were the to- 
tal luminosities measured from the blackbody fitting at each epoch, 
as listed in the final column of Table [7] Table [8] lists the param- 
eters of the best-fitting models for each epoch, including derived 
dust masses, Md, and visual optical depths to. 55- The variation in 
derived dust masses when parameters such as Ri n , Y, L and T are 
varied from their best-fitting value allows us to estimate that the 
uncertainty in the derived dust mass is <50%. Figure [TT] plots the 
emergent SEDs (red solid lines) and the observations, where the 
plotted flux densities have been corrected for foreground extinc- 
tion. 

The observed SEDs were best fitted using diffuse field radi- 
ation temperatures of 7000-8000 K, a little higher than the best- 
fitting hot blackbody temperatures listed in Table [7] This might be 
attributable to the effects of internal dust extinction on the emer- 
gent optical energy distribution in the radiative transfer models. The 
best-fitting models for the different epochs were not homologous, 
in that the R in and R out values did not increase linearly with time 
(Table [8}. Since the derived dust masses increased by a factor of 
four between days 300 and 464, indicating ongoing dust formation, 
the dominant dust-emitting regions could change with time. 

However, we also investigated homologous models for the dif- 
ferent epochs, taking the day 300 value of Y = 4.0 from Tableland 
keeping it the same for subsequent epochs, adjusting only the total 
dust mass to obtain a best fit. These models are plotted as the back 
solid lines in Figure[TT] They produced slightly improved fits to the 
day 406 and 464 24-/xm data-points, although they failed to match 
the observed SED on day 690 (Figure [TTJ. 

Table |8] shows that the smooth model dust masses increased 
from 0.4xl0~ 4 M on day 300 to 1.1 xlO -4 M on day 464 
and 4.4xl0 -4 Mq on day 690. The inner and outer radii of 
the day 300 model correspond to expansion velocities of 2700 
and 10800 km s _1 , respectively, while those for the day 406 
and 464 models correspond to expansion velocities of 2000 and 
8000 km s _1 respectively. The above inner radii velocities are con- 
sistent with line absorption minimum velocities measured in optical 
spectra obtained at these epochs, e.g. |Sahu et al.| p006 ) measured 
Fe II absorption minimum expnasion velocities of ~2000 km s _1 
after day 150. Since absorption line optical depths scale as f ndr, 
where n is the density, then for n oc r~ 2 or steeper density distri- 
butions, line optical depths are strongly weighted to the inner radii, 
where velocities are lowest, as are emission lines, whose emissivi- 
ties typically scale as f Anr 2 n 2 dr. At the earliest epochs however, 
when ejecta densities are much higher, line optical depths of unity 
do not penetrate very deep into the outermost layers of the ejecta, 
where expansion velocities are much higher, e.g. the day 25 spec- 
trum of |Sahu et al.| (2006) showed Ha and Hf3 absorption minima 
at expansion velocities of ^8000 km s _1 (their Fig. 9), consistent 



with the outer radius expansion velocities of our dust models, while 
Ha absorption was even detectable out to —14,500 km s _1 on day 
25. 

7.2 Clumpy dust models 

For the Type II SN 2003gd, |Sugerman et al.|p006] ) demonstrated 
that smooth dust models could underestimate the dust mass by an 
order of magnitude or more compared to models that allow for 
clumping. Ercolano et al.| ( |2007| ) showed that both smooth and 
clumpy dust models could fit the observed SEDs of SN 1987 A 
at late epochs, with clumpy models able to accommodate signifi- 
cantly larger dust masses. We have constructed clumpy models for 
SN 2004et, employing a similar modelling strategy to that used for 
SN 1987 A and SN 2003gd. 

For the clumpy models, we assume that dense homogeneous 
clumps are embedded in a less-dense interclump medium (ICM) 
with an r~ 2 density distribution. The clumps have radius 5xR ou t 
and a volume filling factor, /. We adopted the same S = 1/30 and / 
= 0.01 for all epochs. The density contrast between the clumps and 
the smooth ICM is defined by a = N c i urnp (R in )IN srnoo th(Rin), 
where N c i ump (Rin) and N smooth (Ri n ) are the densities of the 
clumps and the smooth ICM at the inner radius, respectively. We 
set a = 55 for all epochs. For the other parameters, we adopted 
the same source luminosity, temperature, dust composition and size 
distribution as used for the smooth dust models. Similar SEDs to 
those from the smooth dust models were obtained for the clumpy 
models (red dashed curves in Figure[TTJ. Compared to the variable- 
Y smooth dust models, the counterpart clumpy models produced an 
improved fit to the 3.6-/xm photometry on days 300-360 and to the 
24-/xm photometry on days 300-406. 

Compared to the smooth dust distribution models, the 
clumped dust models were able to accommodate two to five times 
larger dust masses without increasing the effective dust optical 
depths in the visible region of the spectrum. The smooth and 
clumpy dust models for day 690 both predict an effective optical 
depth of 1.3 in the V-band, consistent with our estimate from the 
observed light curve of 0.8-1.5 magnitudes of internal extinction 
at this epoch (Section [5^2| . The clumped dust model for 690 listed 
in Table [8] had a total dust mass of 1.5 xlO" 3 M , but we found 
that up to 4x 10 -3 M of dust could be accommodated in clumps 
at that epoch without seriously reducing the goodness of fit to the 
observed SED. 



8 DISCUSSION AND CONCLUSIONS 

We have used a range of new and archival optical and infrared data, 
extending from day 64 to day 2151, to investigate the formation 
of dust in the ejecta of the Type II-P SN 2004et, focusing in par- 
ticular on the mid-IR observations obtained with the Spitzer Space 
Telescope and Gemini-North telescope. 

As discussed by |Lucy et aT](T^89] > for SN 1987A, there are 
three distinct signatures of dust formation in the ejecta of core- 
collapse SNe. These are (1) the appearance of asymmetric blue- 
shifted emission lines caused by dust preferentially extinguishing 
emission from the receding (redshifted) gas; (2) a drop in visual 
brightness due to increased extinction by the newly formed dust; 
accompanied by (3) a mid-infrared excess due to thermal dust emis- 
sion, These three phenomena have all been observed in the case of 
SN 2004et, as found by previous studies ( |Sahu et al.|2006||Kotak| 
|et al.|2009|[Maguire et al.|20T0] ). 
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Figure 11. MOCASSIN dust model fits to the optical-IR (0.3-30 /im) SEDs of SN 2004et between days 300 and 690. Filled circles indicate fluxes observed 
at the epoch of the IRAC observations, whereas open circles indicate optical, NIR and mid-IR photometry which has been interpolated/extrapolated to the 
epochs of the IRAC observations. The filled triangles indicate the Gemini-Michelle N'-band observations, and upper limits to the flux densities are shown as 
downward-pointing arrows. Spitzer IRS spectra obtained at epochs close to the IRAC observations (on days 294, 348 and 480, respectively: green solid lines) 
are in good agreement with the mid-IR photometry. Model fits to the data are indicated by the black solid curves (smooth dust distribution, fixed Y = 4.0), red 
solid curves (smooth dust distribution, varying Y) and the red dashed curves (clumpy dust distribution, varying Y). All models adopted a composition of 20:80 
per cent amorphous carbon: silicates by mass. See text for further details. 



From an analysis of new and recalibrated spectra, we showed 
that between days 259 and 646 the peak of the Ha emission line 
shifted to the blue by 600 km s _1 (Fig. 2). The optical light curve of 
SN 2004et declined more rapidly than that expected from radioac- 
tive deposition, with the onset of ejecta dust formation estimated to 



have occurred ~ 300-400 days after explosion (Fig. 6) The same 
light curves allowed us to estimate that by day 690 the additional 
extinction in the F-band, attributable to newly formed ejecta dust, 
was between 0.8 and 1.5 magnitudes. From day 300 onwards the 
SN 2004et SEDs exhibit clear excess mid-IR emission relative to 
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blackbodies extrapolated from the optical. The day 300, 360 and 
406 SEDs could each be fitted by two blackbody components: (i) 
a hot component attributed to emission from optically thick gas, 
and (ii) a warm component attributed to dust freshly synthesised 
and radioactively heated in the SN ejecta (Fig. 8). While these two- 
component fits provided adequate matches to the observed SEDs 
out to 16 /jm, they did underpredict the 24-/xm fluxes, which re- 
mained relatively constant between days 315 and 709 (Fig. 4). The 
day 294, 348 and 480 Spitzer IRS spectra (Fig. 1 1) also confirm the 
presence of relatively steady excess emission longwards of 20 /im 
at these epochs. By day 464 a cool third component was definitely 
required to fit an increasing excess seen in the 16-/xm photometry. 

For days 300 to 690, our ejecta dust radiative transfer mod- 
elling was able to match the observed SEDs out to 20 /mi (Fig. 11), 
longwards of which the relatively invariant excess emission dis- 
cussed above was present. The observed rise in mid-IR fluxes after 
1000 days coincided with a flattening of the optical and NIR light 
curves at around this time. This, coupled with the fact that the mini- 
mum luminosities estimated from blackbody fitting exceeded those 
expected from radioactive decay from c. day 690 onwards, implies 
that an additional emission source is required after that date. One 
possibility, discussed by Kotak et al. (2009), is that this emission 
source was due to the formation of a cool dense shell of dust as a 
result of ejecta-CSM interaction. Another, discussed by Sugerman 
et al. (in preparation) is that this emission was due to a light echo 
from pre-existing CSM dust. 

For days 300, 360 and 406, the 80% silicate, 20% amorphous 
carbon dust masses derived from our smooth dust models (Table 8) 
agree with the silicate dust masses derived by Kotak et al. ( 2009) 
but our derived dust masses increase faster than their values from 
day 464 onwards and are a factor of 5 larger by day 690. The 
clumped dust models of Kot ak et aL] ([2009) had similar masses 
to their unclumped dust models, while our clumped models have 
dust masses that are factors of 2-4 larger than for our smooth dust 
models (Table 8). The dust mass derived from our day 464 clumped 
model was 5x 10~ 4 M , increasing to 1.5 x 10~ 3 M by day 690. 
This mass of newly formed dust is similar to, or larger than, val- 
ues derived for a number of other recent Type II-P SNe at similar 
epochs, e.g. SN 1987A ([Wooden et al.| 19931 |Ercolano et al.|2007| ), 
SN 2003gd ([Sugerman et al.|2006[|Meikle et al.|2007), SN 2004dj 
JSzalai et al.|2011[|Meikle et al.|201l| >, SN2007it ( [Andrews et"aT] 
|2011| ) and SN2007od ( [Andrews et al.|2010| ). But this mass is still 
a factor of 100 or more smaller than the CCSN ejecta dust masses 
estimated to be needed to account for the large quantities of dust 
observed in some high redshift galaxies ( Kozasa et al. 1989 ; Todini 
|& Ferrara|200T||Nozawa et al.|2003[|Bianchi & Schneider|2007| 
|Dwek et al.|2007l|Dwek & Cherchn eff 201 1) and implies that SNe 
of this type cannot make a major contribution to the dust content of 
galaxies, unless their dust masses continue to grow at later epochs 
than have typically been observed at mid-IR wavelengths. 
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APPENDIX A: PHOTOMETRIC DATA PROCESSING 
Al Mid-IR data processing 

The earliest Spitzer data for the field around SN 2004et were ob- 
tained by the SINGS Legacy team ( Kennicutt et al.|2 003), who pro- 
duced enhanced data products which could be downloaded from the 
SINGS Legacy Data Deliveries section^] of the SSC website. The 
MIPS 24 /xm mosaic of NGC 6946 (providing a pre-explosion im- 
age of the field around SN 2004et at day -75) were downloaded 
from the SINGS fifth and final data delivery. The mosaic was cre- 
ated from multiple Spitzer images obtained in scan-mapping mode 
over two days (75 and 73 days prior to the explosion of SN 2004et) 
and processed with the MIPS Data Analysis Tool version 3.06 
( Gordon et al. 2005), along with additional custom processing by 
the SINGS team. 

However, the enhanced SINGS IRAC data were not used, due 
to the unusually long time span between the two observations of 
NGC 6946 (rsj six months, compared to the usual one or two days) 
which were combined by the SINGS team to construct the final en- 
hanced mosaic. Instead, the standard BCD pipeline data for each 
day (corresponding to the pre-explosion image at day -104 and 

4 http://ssc.spitzer.caltech.edu/legacy/singshistory.html 
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the first mid-IR post-explosion image at day 64) were downloaded 
from the Spitzer archive and processed along with the rest of the 
Spitzer observations as described below. 

The basic calibrated data (BCD) from the Spitzer pipeline 
were combined into final mosaic images using the Spitzer MOPEX 
software package (Makovoz et al. 2006 ), which includes outlying- 
pixel rejection, background matching, and mosaicing with driz- 
zling to increase the sampling of the point-spread function (PSF). 
The IRAC and MIPS data were re-drizzled to plate scales of 
O'.'75/pixel and l'.'5/pixel, respectively (compared to the standard 
Spitzer pipeline products with a spatial sampling of l'.'2/pixel for 
IRAC and 2 / /45/pixel for MIPS), following the pixel sizes adopted 
by the SINGS team for producing their enhanced data products. 
The PUI data were re-drizzled to a plate scale of l'.'2/pixel (com- 
pared to the standard pipeline data with a spatial sampling of 
l'.'8/pixel). The final mosaic images for all three Spitzer instruments 
were calibrated in surface brightness units of MJy sr _1 during the 
BCD pipeline stage. 

For the IRAC and MIPS 24 fim data, PSF-matched difference 
images, which used the pre-explosion SINGS mosaics as the refer- 
ence images, were produced as follows. The final mosaics produced 
by MOPEX were geometrically registered to a common reference 
frame using matching point sources identified within the fields 
of view, with a 2nd-order general fit within the IRAF geotran 
task. In all cases, registration residuals were less than 0.1 pix- 
els RMS in both the x and y dimensions. Once registered, the 
data were PSF-matched and differenced using the DIFIMPHOT 
package (Tomaney & Crotts 1996 ) as implemented and modified 
by |Sugerman et al.| ( 2005| ). Two approaches were taken to PSF- 
matching. In the first, an empirical PSF was built for each image us- 
ing uncrowded point sources combined using the daophot PSF- 
building tasks ( Stetson 1987 ). In the second, these same tasks were 
run on a single theoretical PSF available for each Spitzer IRAC and 
MIPS 24 /am image from the archive, after that model had been ro- 
tated and scaled according to the data's particular plate scale and 
position angle. Images were PSF-matched separately using the em- 
pirical and theoretical PSFs, photometrically scaled by the median 
brightness of a number of matching point sources, and then sub- 
tracted to yield the final difference images. In general, the differ- 
ence images made using the theoretical PSFs were of higher quality 
{i.e., smaller subtraction residuals and less background noise) since 
the shapes of Spitzer' 's PSFs have varied quite little during its mis- 
sion. In practice, these techniques allow the reliable detection and 
measurement of changes in point sources that would be considered 
below the l-a level in direct, undifferenced images ( Sugerman et al. 
12002] ). 

For the photometry from the IRAC and MIPS difference im- 
ages, uncertainties were measured with a custom implementation of 
an optimal photometry code written for the original version of DI- 
FIMPHOT, which includes the full noise model of the daophot- 
allstar routine as well as correct noise contributions from both 
input and reference images and the original sky values in each im- 
age prior to pipeline calibration. 

The Gemini Michelle data were downloaded from the 
Gemini Science Archive and processed with the Gemini 
IRAF midir tasks and further cleaning procedures. The fluxes, or 
upper limits, in counts measured from the final average-combined 
images were converted to F u flux density units by multiplying by 
flux conversion factors derived from aperture photometry of stan- 
dard stars. 



A2 Optical and near-IR data processing 

Gemini/GMOS-N On days 317, 404, and 664, 60s images were 
taken with GMOS-N in the g f ,r f , and i' broad-band filters. Longer 
exposures of 2 x 600 s in the same filters were taken at day 1412 
when the SN was expected to have faded. The images were reduced 
using the Gemini IRAF package. Pipeline processed bias and flat 
field images were obtained from the Gemini Science Archive. Ob- 
ject images were trimmed, corrected for overscan and bias, and flat- 
fielded using the gs reduce task. Finally, gmosaic was used to 
mosaic the three GMOS CCDs into a single image. PSF-fitted pho- 
tometry was performed on the SN and a sample of standard stars 
from the photometric V, R, and I sequence of Pozzo et al. ( 2006]) to 
establish the nightly zeropoint. The GMOS Sloan Digital Sky Sur- 
vey (SDSS) magnitudes were transformed into Johnson V, R, and 
/ magnitudes using the transformation equations given by Welch 
|et al.| ( |2007| >. These linear transformations were derived from the 
photometric V, R, and I sequence presented by |Pozzo et al. (2006). 
Gemini/NIRI In order to ensure a consistent sky background be- 
tween exposures, individual exposures of 30 s in each of the broad- 
band JHK filters were taken, with the total number being dictated 
by the anticipated decline of the SN NIR light. A 5 " dither pattern 
was employed to ensure efficient removal of point sources while 
making sky images. Data reduction for each night was performed 
using the standard NIRI routines within the Gemini IRAF pack- 
age, nprepare and nif lat were used to derive the normalized 
flat field and the bad pixel mask while nisky was used to cre- 
ate the final sky image, ni reduce was used to subtract this sky 
image from and apply the flat field correction to the processed ob- 
ject images. Finally, the individual images in each filter were coad- 
ded using the GEMTOOLS routine imcoadd. PSF-fitted photom- 
etry was performed on SN 2004et and three standard stars present 
within the field. The JHK magnitudes of the standard stars are 
contained within the 2MASS catalogue and were used to derive 
the nightly zeropoint in each filter. The photometric uncertainty is 
dominated by the standard deviation of the zeropoint derived from 
the three standard stars. 

Tenagra and Steward/Bok These early optical and NIR data 
were reduced, calibrated and measured in a similar manner to the 
GMOS-N and NIRI observations described above, using standard 
routines within IRAF. The data were bias and dark subtracted, flat- 
fielded, and the multiple exposures were combined to form the final 
image. 

HST WFPC2 and NICMOS2 For the WFPC2 observations, imag- 
ing with the broadband F606W and F814W filters was used to 
measure the V- and /-band flux densities. For the NICMOS2 ob- 
servations, imaging was carried out with the F110W, F160W, and 
F205W filters, most closely representing the standard JHK-bmds. 
Small-scale dithering was employed for both instruments to im- 
prove S/N, remove cosmic rays and to improve the pixel- scale of 
the final images by drizzle techniques. The WFPC2 observations 
used a 4-point dither x 400 s exposure for each band. The NIC- 
MOS2 observations used a 5-point dither x 128 s in the F110W 
and F205W filters and a 4-point dither x 128 s for F160W. 

The data were reduced and calibrated using the IRAF exter- 
nal package stsdas (version 3.8) and included the removal of 
cosmic rays and other artifacts, as well as linearity corrections. 
High-resolution images were created using the s t s da s/dr i z z 1 e 
package and additional distortion correction and alignment was 
performed using background stars. The point spread function of 
these reference background stars were fit by Gaussian profiles 
to obtain accurate positions, and then instrumental distortions 
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were corrected with the IRAF tasks xyxymatch, geomap, and 
geotran. The resultant pixel scale is ^0.02" pixel -1 in the 
WFPC2 images and ~0.04" pixel -1 in the NICMOS2 images. The 
FWHM of the PSFs was ~3 pixel. Flux measurements were per- 
formed using the IRAF da o phot tasks. For the WFPC2 data, 
HST magnitudes were converted to the Johnson-Cousins system us- 
ing the transforms of Dolphin ( 2000 2009 ), which include charge 
transfer efficiency (CTE) corrections. 

Whilst the broad band NICMOS filters, F110W, F160W and 
F205W, are roughly equivalent to the J, H and K filters respec- 
tively, colour transformations are not well constrained for late-time 
SN spectra at NIR wavelengths, and consequently magnitudes were 
not converted to the standard JHK photometric system. The late- 
time ffiT-NICMOS magnitudes listed in Tableware for the HST 
filters in the Vegamag system, which uses an estimate of the flux 
density of Vega, from synthetic spectra integrated over the NIC- 
MOS bandpasses, as a photometric zeropoint. For each filter and 
epoch, the measured count rate (CR, in units of DNs -1 ) at the 
position of the SN was converted to flux by multiplication with 
the PHOTFNU (JysDN -1 ) conversion factor given in the fits 
header, where PHOTFNU is the bandpass-averaged flux density 
for a source that would produce a count rate of 1 DN -1 . An approx- 
imate Vega-normalised magnitude was then calculated from the fol- 
lowing equation (NICMOS data handbook; [Thatte et al.|2009| >: 

m = ZPtVeg^^.Slog^iPHOTFNUxCRxiFviVega))- 1 

where (F u (Vega)) is the bandpass averaged flux density (in Jy) 
for the NICMOS filters using a model reference spectrum of Vega[^] 
(Bohl in|2007t and ZP(Vega) is the magnitude of Vega, which is 
defined to be 0.00 mag under the California Institute of Technology 
(CIT) infrared photometry scale. 

WIYN/WHIRC The WHIRC ( |Meixner et "all |20T0] ) has a 
2048x2048 HgCdTe VIRGO detector with a pixel scale of - 07 1 
pixel -1 . Sky conditions during the observation were fair, with a 
seeing of ~ 0"8. To minimise the effects of high background lev- 
els and pixel-to-pixel variations on the array, dithering techniques 
were employed whereby the source was offset in each frame of a 
series of exposures. The offset images were used for sky level cor- 
rections. Data reduction was carried out using standard IRAF tasks. 
The array linearity correction was performed using the WHIRC 
task wprep and a distortion correction was applied using files 
downloaded from the WIYN- WHIRC web pagqj A selection of 
2MASS stars close to SN 2004et were used for final flux cali- 
bration, and flux measurements were performed using the IRAF 
daophot package. 

The WIYN if -band observation of SN 2004et on day 1803 
was the last of our NIR observations of the SN, some 5 years af- 
ter explosion. The SN was not clearly detected due to the con- 
tribution from neighbouring stars that were resolved in the high- 
resolution ffiT-NICMOS observations (on days 1019 and 1215; 
see the following section). Consequently, an upper limit to the mag- 
nitude was derived. The magnitude in an aperture of radius 0"6 
(~ 17 pc, for the adopted distance of 5.9 Mpc) was measured using 
the IRAF daophot tasks. Estimated contributions from 3 neigh- 
bour stars resolved in the high-resolution NICMOS F160W filter 



5 Taken from the NIC2 table of Photometric Keywords and Vegamag 
Zeropoints at 

http://www.stsci.edu/hst/nicmos/performance/photometry/postncsJceywords.html 
° http://www.noao.edu/wiyn/ 



(~ H band) were measured from the NICMOS data and subtracted 
from the WIYN magnitude to provide the final upper limit pre- 
sented in Table |6] 



A3 Late-time high resolution HST images 

The late-time high resolution HST images reveal that the single 
point source seen at the SN position in the Spitzer images is actually 
a complex field comprised of at least 3 sources. Figure [AT] shows 
example HST WFPC2 F606W (« V) and NICMOS F110W 
J) images from July 2007 (day 1019) compared to the Spitzer 
IRAC 3.6-fim image from August 2007 (day 1054). Each field of 
view is centred on the position of the SN and a source is located 
at this position indicated by the cross-hairs. The star to the east 
("star 2") was detected in both WFPC2 filters and all three NIC- 
MOS filters at days 1019 and 1215. The third star to the south of 
the SN position ("star 3") was detected in all three NICMOS bands 
but was not detected in the WFPC2 filters at either epoch. 

The magnitude of star 2 was measured from the WFPC2 im- 
ages and transformed to V and I c band magnitudes as described 
previously for the SN. 

The V and I c magnitudes of star 2 were de-reddened using 
E(B -V) = 0.41 mag dZwitter et al.|2004| ) and adopting the extinc- 
tion law of |Cardelli et al.| ( 1989| > withitV = 3.1, corresponding 
to Ay — 1.27 ± 0.22 mag. The intrinsic (V — I) c colour was 
then used to estimate an R c band magnitude using a table of in- 
trinsic colours as a function of spectral type compiled by t he Space 
Telescope Science Institute (STSclFJbased on the work of Fitzger- 
ald (1970) and Ducati et al. (2001). The derived colour index of 
(V — I) c — 0.8 indicated star 2 was probably a K1.0 star with 
an intrinsic R c magnitude of ^22.5. The equivalent reddened R c 
magnitude of 23.53 =t 0.50 was in reasonable agreement (within 
the errors) with a magnitude estimated from fitting two PSFs to 
the blended 'SN plus star 2' in the i? c -band GMOS-N image at day 
1412. This PSF-fitted measurement gave an R c band magnitude for 
star 2 of 23.63 ± 0.19. A B band magnitude for star 2 of 25.46 ± 
0.50 was estimated in the same way. 

The estimated reddened B and R c magnitudes of star 2, to- 
gether with the measured V and I c magnitudes, were subtracted 
from the optical magnitudes of SN 2004et obtained from the 
Subaru-FOCUS and GMOS-N observations from day 646 onwards. 
During this time the SN had faded substantially such that the neigh- 
bouring star 2 made a significant contribution to the brightness 
measured at those epochs. At day 646, the estimated contribution 
from star 2 was ~ 8 % in the V band and ~ 10 % in the I c band, 
whilst by day 1412, the contribution was almost half that of the total 
flux in both bands. For epochs earlier than day 500, the brightness 
contribution from star 2 was < 1 %. 

The late-time Gemini-NIRI photometry at day 652 has not 
been corrected for contamination by stars 2 and 3. As previously 
explained, the NICMOS magnitudes are listed in the HST Vegamag 
system, and were not converted to the standard JHK photometric 
system since colour transformations are not well constrained for 
late-time SN spectra at NIR wavelengths. It is possible that these 
neighbouring stars made a small but significant (~ 10%) contri- 
bution to the brightness measured in the NIRI observations at day 
652, so the uncertainties on the day 690 interpolated flux densities 
in Table [6] have been increased to reflect this, by adding a further 
10 % error in quadrature to the original flux uncertainty. 



http://www.stsci.edu/ inr/intrins.html 



24 Fabbri et al. 









i i # i | i i | , , , | , 


1 1 




_ 1 1 1 


*j I i i i | i i i | i 


1 1 _ 












4 








4 








4 














1 








- "•'>.' 








■ 




















2 








2 




* 




2 












1 






o 

(D 

m 




































u 

D 





























1 












































< 










-2 
























-2 














-2 
















































-4 




(a) WFPC2 F606W 

i i i 1 i i i 1 i i i 1 i 


1 , " 


-4 


I (b) NICMOS F1 1 OW 

~ i I i i i I i i i I i i i I i 


1 , 1 


-4 


■ (c 

"i 1 


) IRAC 3.6 


, , , 1 




1 








-4-2 2 


4 




-4 


-2 2 


4 




-4 




2 




2 




4 





Aa (arcsec) Aa (arcsec) Aa (arcsec) 



Figure Al. Late-time high-resolution HST images reveal a complex field compared with Spitzer IRAC data. Panel (a) shows the WFPC2 F606W image at 
day 1019 with a 5" x5 /; FOV centred on the SN position. A source assumed to be the SN is located at the centre of the image (indicated by the cross-hairs), 
with a close companion to the east ("star 2") which is detected in both WFPC2 filters. Panel (b) shows the equivalent field in the NICMOS F110W filter at the 
same epoch. Star 2 is detected in all 3 NICMOS filters together with an additional red object just south of the SN position. For comparison, panel (c) shows 
the Spitzer IRAC 3.6-/zm image at day 1054. 



